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ABSTRACT 

Elastothermodynamic  Damping  in  Composite  and  Cracked  Media. 

(December  1995) 

Joseph  Edward  Bishop,  B.S.,  Texas  A&M  University; 

M.S.,  Texas  A&M  University 
Chair  of  Advisory  Committee:  Dr.  Vikram  K.  Kinra 

When  a  composite  material  is  subjected  to  a  stress  field,  different  regions 
undergo  different  temperature  fluctuations  due  to  the  well-known 
thermoelastic  effect.  As  a  result  irreversible  heat  conduction  occurs,  and 
entropy  is  produced  which  is  manifested  as  a  conversion  of  mechanical  energy 
into  heat.  Moreover,  the  changes  in  temperature  produce  a  thermal  strain  that 
is  out  of  phase  with  the  stress,  thus  converting  mechanical  energy  into  heat, 
i.e.  work  is  lost.  We  define  this  process  as  elastothermodynamic  damping. 

Herein,  using  the  linear  one-way  coupled  theory  of  elastothermodynamic 
relaxation,  the  elastothermodynamic  damping  of  composite  and  cracked  media 
is  examined.  Two  equivalent  descriptions  of  elastothermodynamic  damping 
are  first  established:  (1)  the  mechanical  description  and  (2)  the  entropic 
description.  An  integral-transform  technique  is  then  developed  to  solve  for 
the  temperature  field  of  a  general  composite  medium  with  arbitrary  heat 
generation.  With  this  solution,  a  general  expression  for  the  damping  of  the 
composite  material  is  derived.  This  general  solution  for  the 
elastothermodynamic  damping  is  then  specialized  to  a  solution  for  a 
composite  material  consisting  of  A/-isotropic  layers  in  a  rectangular, 
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cylindrical,  and  spherical  coordinate  system  (i.e.  an  A^-layer  slab,  cylinder,  and 
sphere,  respectively),  subjected  to  any  stress  field  so  long  as  the  resulting  heat 
conduction  can  be  described  by  a  single  spatial  coordinate  orthogonal  to  the 
layering.  With  this  specialized  result  the  elastothermodynamic  damping  is 
examined  in  laminated,  fiber-reinforced,  and  particulate  composites.  Finally, 
the  elastothermodynamic  damping  of  cracked  media  is  examined. 
Specifically,  an  approximate  analysis  is  given  for  the  entropy  produced  and 
work  lost  in  the  neighborhood  of  a  Griffith  crack  subjected  to  a  time-harmonic 
loading  in  Modes  1, 11,  and  III.  In  all  three  Modes  the  temperature  at  the  crack 
tip  remains  bounded.  In  Mode  I  the  entropy  produced  (per  unit  volume  per 
cycle)  is  finite  at  the  crack  tip,  whereas  the  work  lost  (per  unit  volume  per 
cycle)  goes  to  infinity  as  l/Vr .  Conversely,  in  Mode  II  the  entropy  produced 
goes  to  infinity  as  1/r  as  the  crack  tip  is  approached,  whereas  the  work  lost  is 
finite.  In  Mode  III  the  thermoelastic  effect  disappears  altogether  and, 
therefore,  both  the  entropy  produced  and  the  work  lost  are  zero  throughout  the 
plate. 
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1.  INTRODUCTION 

Composite  materials  are  increasingly  the  material  of  choice  for  a  wide 
variety  of  engineering  applications.  Whereas  a  considerable  amoimt  of 
research  has  been  performed  concerning  the  stiffiiess,  strength,  fracture 
toughness,  fatigue  life,  coefficient  of  thermal  expansion,  etc.,  of  composite 
materials,  their  damping  properties  have  been  investigated  to  a  much  lesser 
extent.  However,  damping  is  a  critically  important  material  property  from  the 
viewpoint  of  vibration  suppression  in  a  variety  of  applications:  submarines, 
aircraft  structures,  automobile  bodies,  office  machines,  off-shore  structures, 
and  in  particular,  large  and  flexible  aerospace  structures,  e.g.  a  space  station. 
A  significant  fraction  of  the  body  of  literature  on  damping  concerns  itself  with 
viscoelastic  damping  in  polymeric  materials  and  polymer-matrix  composites 
(PMCs).  Even  though  these  materials  exhibit  high  damping,  there  are  many 
situations  where  they  caimot  be  used,  for  example,  elevated  temperature 
applications  such  as  the  high-speed  civil  transport,  turbine  engines,  etc. 
Moreover,  they  cannot  be  used  in  space  structures  due  to  the  problem  of 
degassing  and  an  attendant  degradation  of  mechanical  properties.  Therefore, 
there  is  a  need  for  developing  new  materials  with  high  intrinsic  material 
damping  without  significantly  sacrificing  other  mechanical  properties.  As  a 
result,  there  has  recently  been  a  considerable  activity  in  designing  and 
characterizing  high-damping  metal-matrix  composites  (MMCs)  (for  example. 

This  dissertation  follows  the  style  and  format  of  the  Journal  of  Applied 
Mechanics. 
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DiCarlo,  1976;  DiCarlo  and  Maisel,  1979;  Bhagat,  et  al,  1989;  Kinra,  et  al, 
1991;  Kinra  and  Wren,  1992;  Wren  and  Kinra,  1992;  Wong  and  Holcomb, 
1992). 

For  a  given  material  there  are  many  mechanisms  contributing  to  the  total 
damping  of  an  MMC,  e.g.  point  defect  relaxation,  dislocation  motion,  grain 
boundary  sliding,  coulomb  friction  at  the  inclusion-matrix  interface, 
magnetoelastic  effects,  and  elastothermodynamic  relaxation  (Nowick  and 
Berry,  1972;  De  Batist,  1972;  Kinra  and  Wolfenden,  1992).  In  a  recent  book- 
chapter  Wolfenden  and  Wolla  (1992)  have  done  an  excellent  review  of 
damping  in  MMCs  .  From  the  results  of  these  investigations,  it  has  been 
found  that  most  of  the  damping  mechanisms  offer  a  significant  damping  only 
over  a  small  range  of  frequency,  temperature,  and  stress  (or  strain). 
Elastothermodynamic!  damping  (also  known  as  thermoelastic  damping) 
offers  an  attractive  alternative.  In  this  case  the  damping  is  due  to  thermal 
currents  within  the  material  and  is  not  accompanied  by  any  significant  loss  of 
stiffness  or  strength.  To  first  order  it  is  linearly  dependent  on  the  absolute 
temperature,  and  it  is  independent  of  the  stress  level  (i.e.  damping  is  linear). 
As  we  shall  see,  analogous  to  other  sources  of  damping,  elastothermodynamic 
damping  also  suffers  from  the  limitation  that  it  is  significant  only  over  a  finite 


*  We  propose  that  elastothermodynamic  is  a  more  appropriate  qualifier  for  the  word  damping 
than  thermoelastic.  Historically  (Boley  and  Weiner,  1960),  the  word  thermoelastic  has  been 
used  to  connote  a  situation  where  the  cause  is  an  externally  applied  temperature  field 
{thermo),  and  the  effect  is  an  elastic  stress  field  {elastic),  hence  thermoelastic.  In  the  present 
context  of  damping,  there  is  no  externally  applied  temperature  field.  Instead,  the  cause  is  an 
externally  applied  traction  field  and  the  accompanying  stress  field  (elasto),  and  the  effect  is  a 
temperature  field  {thermo).  Moreover,  the  dissipation  of  mechanical  energy  occurs  only  if  the 
stress  field  is  time  dependent  {dynamic),  hence  elastothermodynamic 
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range  of  frequency.  Finally,  it  is  envisioned  that  the  elastothermodynamic 
damping  of  composites  can  be  tailored  in  a  manner  analogous  to  that  used 
currently  for  tailoring  stiffness,  strength,  coefficient  of  thermal  expansion,  etc. 
Of  course  the  tailoring  of  the  former  cannot  be  done  independently  of  the 
tailoring  of  the  latter. 

Whenever  a  material  is  stressed  adiabatically,  there  is  an  accompanying 
change  in  temperature,  however  small.  This  phenomenon  is  known  as  the 
thermoelastic  effect  and  can  easily  be  deduced  from  the  Maxwell's  relations 
(e.g.  Zemansky  and  Dittman,  1981).  For  an  anisotropic  material  this  effect  is 
given  by 


_ _ ’J’^USL 


(1) 


where  T  is  the  absolute  temperature,  is  the  stress  tensor,  a,„  is  the  linear 
coefficient  of  thermal  expansion  tensor,  C  is  the  specific  heat  per  unit  volume, 
and  the  subscript  5  denotes  an  isentropic  process.  Since  the  temperature  and 
stress  fields  are  coupled  via  Eq.  (1),  inhomogeneities  in  stress  and/or  material 
properties  will  result  in  inhomogeneities  in  temperature.  Heat  will  conduct 
from  the  high-temperature  regions  to  the  low-temperature  regions,  and  as  a 
consequence  of  the  resulting  irreversible  heat  conduction,  entropy  will  be 
produced.  In  a  cyclic  process  such  as  a  vibrating  structure,  this  entropy 
production  will  be  manifested  as  a  conversion  of  external  work  into  internal 
heat  energy.  We  define  this  process  as  elastothermodynamic  damping.  This 
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increase  in  the  internal  heat  energy  over  one  cycle  of  vibration  will  also  be 
referred  to  as  the  work  lost. 

The  measurement  of  material  damping  is  a  notoriously  difficult  problem. 
Any  irreversible  phenomenon  will  increase  the  measured  damping,  e.g.  air 
damping,  inelastic  clamping  devices,  etc.  Furthermore,  as  mentioned 
previously,  there  are  many  mechanisms  of  material  damping,  and  isolating 
any  one  mechanism  for  study  is  experimentally  demanding.  The  first 
experimental  research  into  elastothermodynamic  damping  was  performed  by 
Zener  in  1938(b)  using  a  copper  beam.  He  observed  that  the  measured 
damping  agreed  with  the  theory  only  with  respect  to  the  shape  of  the  damping 
versus  frequency  curve  and  with  respect  to  the  fi-equency  of  peak  damping.  A 
precise  verification  of  elastothermodynamic  damping  as  a  function  of 
fi-equency  was  performed  by  Berry  in  1955  using  beams  made  of  a-brass.  The 
comparison  between  theory  and  experiment  was  excellent.  Most  recently, 
elastothermodynamic  damping  has  been  measured  as  a  function  of  frequency 
in  6061-T6  aluminum  beams  by  Wren  and  Kinra  (1988)  and  Bishop  and  Kinra 
(1992)  with  an  excellent  comparison  between  theory  and  experiment. 

Elastothermodynamic  damping  has  been  studied  analytically  by  several 
researchers  over  the  past  sixty  years.  Following  Zener's  seminal  work  (1937, 
1938a),  many  authors  have  studied  this  damping  mechanism  in  homogeneous 
structures  using  a  variety  of  solution  techniques  (Biot,  1956;  Liicke,  1956; 
Deresiewicz,  1957;  Alblas,  1961,  1981;  Chadwick,  1962;  Goodman,  et  al., 
1962;  Tasi,  1963;  Tasi  and  Herrmann,  1964;  Gillis,  1968;  Shieh,  1971,  1975, 
1979;  Lee,  1985;  Landau  and  Lifshitz,  1986).  Budiansky  and  Sumner  (1983) 
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studied  bvilk  attenuation  of  composite  materials  using  the  self-consistent 
scheme.  Armstrong  (1984)  analyzed  the  attenuation  of  a  longitudinal  wave 
passing  through  either  a  periodic  or  a  random  sequence  of  slabs.  Eshelby 
(1949)  was  the  first  to  study  the  elastothermodynamic  damping  of  an  internal 
defect,  namely  a  quasistatic  edge  dislocation.  Weiner  (1958)  extended 
Eshelby 's  work  to  high-speed  edge  dislocations.  Savage  (1966)  studied  the 
elastothermodynamic  attenuation  of  a  medium  containing  elliptical  cavities. 

The  objective  of  this  dissertation  is  to  analytically  study  the 
elastothermodynamic^  damping  in  composite  and  cracked  media  subjected  to 
time-harmonic  tractions.  The  linear  one-way  coupled  theory  of 
elastothermodynamic  damping  is  used  for  all  analysis.  In  Section  2  the  linear 
one-way  coupled  theory  of  elastothermodynamic  damping  is  presented.  In 
Section  3  the  elastothermodynamic  damping  of  composite  media  is  examined. 
In  Section  3.1  two  descriptions  of  elastothermodynamic  damping  are 
established:  (1)  the  mechanical  description  and  (2)  the  entropic  description. 
In  Section  3.2  the  solution  of  the  inhomogeneous  heat-conduction  problem 
(the  heat-conduction  problem  including  heat  generation)  by  the  integral- 
transform  technique,  given  by  Ozisik  (1980,  pp.  594-610)  for  a  composite 
medium  consisting  of  isotropic  subregions  with  thermally  imperfect 
interfaces,  is  extended  to  a  composite  medium  consisting  of  several 
subregions,  each  of  which  may  be  anisotropic.  Using  the  results  of 
Sections  3.1  and  3.2,  a  general  expression  for  the  damping  is  derived  for  the 

^Henceforth,  the  term  damping,  when  used  without  qualification,  will  imply 
elastothermodynamic  damping. 
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composite  material  in  Section  3.3.  In  Section  3.4  this  general  solution  for  the 
elastothermodynamic  damping  is  then  specialized  to  a  solution  for  a 
composite  material  consisting  of  iV-isotropic  layers  in  a  rectangular, 
cylindrical,  and  spherical  coordinate  system  (i.e.  an  A^-layer  slab,  cylinder,  and 
sphere,  respectively),  subjected  to  any  stress  field  so  long  as  the  resulting  heat 
conduction  can  be  described  by  a  single  spatial  coordinate  orthogonal  to  the 
layering.  Many  problems  of  current  engineering  interest  can  be  modeled  as 
special  cases  of  the  iV-layer  solution.  (1)  Cartesian  coordinate  system: 
laminated  composites  in  bending  and  extension;  aligned  platelet  composite;  or 
a  periodic  array  of  labs.  (2)  Cylindrical  coordinate  system:  a  unidirectional 
fiber-reinforced  composite  where  each  fiber  may  consist  of  several  concentric 
cylinders.  (3)  Spherical  coordinated  system:  a  particulate  composite  where 
each  particle  may  be  solid  or  hollow,  or  consist  of  several  concentric  spheres, 
for  example,  an  aluminum  matrix  filled  with  hollow  AI2O3  microspheres. 
Specific  results  are  given  for  the  following  problems.  (I)  A  periodic  array  of 
isotropic  slabs  (rectangular  symmetry)  with  a  two-layer  representative  volume 
element  and  a  thermally  imperfect  interface  subjected  to  the  following 
mechanical  stress  states:  (a)  a  time-harmonic  uniform  stress  perpendicular  to 
the  layering,  and  (b)  a  time-harmonic  uniform  strain  parallel  to  the  layering. 
(II)  A  symmetric  three-layer  plate  in  pure  bending  with  isotropic  layers  and 
with  thermally  perfect  interfaces.  (Ill)  A  hollow  two-layer  composite  sphere 
in  which  each  layer  is  isotropic  and  has  a  thermally  imperfect  interface, 
subjected  to  a  uniform  time-harmonic  radial  load.  In  Section  4  the 
elastothermodynamic  damping  of  cracked  media  is  examined.  Specifically,  an 
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approximate  analysis  is  given  for  the  temperature  field,  the  work  lost,  and  the 
entropy  produced  in  the  vicinity  of  a  Griffith  crack  in  Modes  I,  II,  and  III. 
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2.  THE  LINEAR  ONE-WAY  COUPLED  THEORY  OF 
ELASTOTHERMODYNAMIC  DAMPING 


Consider  a  solid  subjected  to  cyclic  forced  vibrations.  If  the  solid  is 
inelastic,  a  portion  of  the  external  work  performed  by  the  applied  tractions 
will  be  converted  into  heat  and  other  nonrecoverable  forms  of  stored  energy. 
Without  reference  to  any  particular  physical  mechanism  of  inelasticity 
(damping),  the  work  lost  per  unit  volume  during  a  cycle  of  vibration  Aw  is 
given  by 

Aw  =  |g,„c/s,„  (2) 

where  is  the  strain  tensor  (Nowick  and  Berry,  1972,  pp.  13-14).  The  work 
lost  per  cycle  AW  over  the  entire  material  occupying  a  region  R  is  given  by 

AW  =  Aw  dR  =  jdRj  (y,„dz,„  (3) 

This  work  lost  may  be  interpreted  as  the  work  required  per  cycle  to  maintain 
the  solid  in  steady  vibration.  A  normalized  measure  of  the  external  work  lost 
is  given  by  the  specific  damping  capacity,  which  is  defined  as 


W 


(4) 


where  W  is  the  maximum  elastic  energy  during  the  cycle. 

Now  consider  a  particular  mechanism  of  anelasticity,  namely, 
elastothermodynamic  damping  (The  words  "inelastic"  and  "anelastic"  are  used 
in  the  sense  defined  by  Nowick  and  Berry,  (1972,  pp.  2-3)).  For  simplicity 
only  a  solid  body  with  adiabatic  outer  boundaries  is  considered.  The  rate  of 
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entropy  production  per  unit  volume,  Sp ,  due  to  irreversible  heat  conduction  is 
given  by 


where  q,  is  the  heat  flux  vector  and  (  )„  denotes  partial  differentiation  with 
respect  to  the  components  of  the  position  vector  x,  (de  Groot  and  Ma2ur, 
1984).  Herein,  it  is  assumed  that  the  heat  flux  vector  satisfies  the  Fourier  law 
of  heat  conduction 

<li  =  -KT,m  (6) 

where  k,^  is  the  thermal  conductivity  tensor  (Ozisik,  pp.  611-629).  The 
entropy  produced  per  unit  volume  As  during  a  cycle  of  vibration  is  given  by 

lxs  =  ^ipdt  (7) 

An  additional  entropy  is  produced  if  the  medium  has  any  thermally 

imperfect  interfaces  at  which  the  temperature  is  discontinuous.  This  entropy 
production  will  be  described  in  Section  3. 1.3 .3.  The  entropy  produced  per 
cycle  A5  over  the  entire  body  occupying  a  region  R  is  thus 

A5=  f  +  (8) 

JR 

The  temperature  field  is  governed  by  the  two-way  coupled  heat 
conduction  equation 

(Nowinski,  1978,  Ch.  6).  The  term  on  the  right-hand  side  is  due  to  the 
thermoelastic  effect,  Eq.  (1).  The  qualifier  "two-way  coupled"  implies  that  the 
stress  field  affects  the  temperature  field,  and  conversely  the  temperature  field 
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affects  the  stress  field.  For  most  crystalline  solids  the  relative  temperature 
change  resulting  from  the  thermoelastic  effect  is  very  small.  For  example,  if 
an  imstressed  rod  of  aluminum  at  an  equilibrium  temperature  of  7^  =  300 AT  is 
subjected  to  a  uniaxial  tensile  stress  of  100  MPa,  the  decrease  in  temperature 
is  0.3  K,  i.e.  |Ar/2^|  =  10'^.  It  is  well  known  that  the  stresses  produced  by 
these  small  changes  in  temperature  are  negligibly  small  compared  to  the 
applied  stresses  (Boley  and  Weiner,  1960).  Therefore,  for  our  purpose  of 
calculating  the  small  changes  in  temperature  resulting  jfrom  the  thermoelastic 
eflfect,  we  assume  that  the  stress  field  is  independent  of  the  temperature  field. 
Eq.  (9)  now  becomes  the  linear  one-way  coupled  heat-conduction  equation. 
Moreover,  without  introducing  an  appreciable  error,  the  algebraic  effort  can  be 
greatly  reduced  by  replacing  T  by  2^  in  the  last  term  (Nowinski,  1978, 
Ch.  6.1);  then  Eq.  (9)  reduces  to 

(10) 

Similarly,  Tis  replaced  with  7^  in  Eq.  (5);  then,  using  Eq.  (6), 

(11) 

These  equations  are  collectively  referred  to  as  the  linear  one-way  coupled 
theory  of  elastothermodynamic  damping.  A  simple  example  illustrating  the 
physics  of  both  the  linear  and  nonlinear  theories  of  elastothermodynamic 
damping  is  given  in  the  Appendix. 


11 


3.  ELASTOTHERMODYNAMIC  DAMPING  IN 
COMPOSITE  MATERIALS 

In  this  Section  the  elastothermodynamic  damping  of  composite  materials 
loaded  time-harmonically  is  examined.  First,  the  equivalence  of  the 
mechanical  and  entropic  descriptions  of  elastothermodynamic  damping  for  a 
general  composite  medium  is  established. 

3.1  EQUIVALENCY  OF  THE  MECHANICAL  AND  ENTROPIC 
DESCRIPTIONS  OF  ELASTOTHERMODYNAMIC  DAMPING  FOR  A 
GENERAL  COMPOSITE  MEDIUM 

In  1938  Zener  showed  that  within  the  linear  one-way  coupled  theory  of 
elastothermodynamic  damping,  AW=T„AS,  where  T^AS  is  simply  the 
increment  of  heat  per  cycle  AQ  in  the  medium.  Thus,  the 
elastothermodynamic  damping,  characterized  by  'F,  may  be  obtained  by 
calculating  either  the  total  external  work  lost  per  cycle  AW  or  the  heat 
increment  per  cycle  T^AS.  The  first  approach  will  be  referred  to  as  the 
mechanical  description  and  the  second  approach  as  the  entropic  description. 
However,  Zener's  analysis  is  confined  to  isotropic  homogeneous  bodies^. 

The  objective  of  this  section  is  to  extend  Zener's  results  to  a  composite 

material  made  up  of  several  subregions,  each  of  which  may  be  anisotropic, 

^Zener's  analysis  is  slightly  more  general  than  what  is  mentioned  here.  In  his  derivation  he 
uses  the  fact  that  for  an  isotropic  material,  ki^  =  and  =  a5;„,  which  is  also  true 
for  a  material  with  cubic  symmetry. 
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and  whose  interfaces  may  be  thermally  imperfect.  In  the  past  Kinra  and 
Milligan  (1994)  have  assumed  that  Aw=  7^ As,  i.e,  the  local  work  lost  at  a 
point  per  cycle  is  equal  to  the  local  entropy  produced  per  cycle  at  that  point 
times  the  equilibrium  temperature.  The  second  objective  of  this  section  is  to 
demonstrate  that  in  general  this  statement  is  not  true.  Since  Zener's  derivation 
was  very  terse  and  extremely  difficult  to  follow,  his  result  is  first  rederived 
using  a  slightly  different  approach  (and  hopefully  easier  to  follow). 
Moreover,  the  framework  developed  in  this  rederivation  will  allow  for  the 
construction  of  a  coimterexample  to  the  statement  Aw  =  T^As. 

3.1.1  Proof  that  AW  =  T^AS  for  a  Homogeneous  Isotropic  Solid 

The  fact  that  for  an  isotropic,  homogeneous,  thermoelastic  solid  with 
adiabatic  boimdaries,  subjected  to  time-harmonic  tractions,  the  total  work  lost 
per  cycle  throughout  the  body,  AW,  is  equal  to  the  total  heat  increment  per 
cycle  throughout  the  body,  7^  AS ,  will  now  be  proven. 

3.1.1.1  Temperature 

For  an  isotropic  material,  k,„  =  where  k  is  the  thermal  conductivity, 
is  the  Kronecker  delta,  and  a,„  =  a5,„,  where  a  is  the  coefficient  of  linear 
thermal  expansion,  Eq.  (10)  reduces  to 

(kT,.),.-C^=Za^  (12) 

ot  ot 

The  stress  field  may  be  separated  into  its  spatial  and  temporal  parts. 


(13) 
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where  x  is  a  position  vector,  co  is  the  frequency  of  vibration.  Herein,  the 
spatial  part  of  the  stress  tensor,  C/„(x),  is  taken  to  be  real  valued.  Since 
Eq.  (12)  is  a  linear  partial  differential  equation,  let 

r(x,o  =  r;+F(x,(D)e'“'  (14) 

where  T*  =(l+z)2^,  and  V  is  the  complex  fluctuation  of  temperature  about 
T* .  As  usual  a  physically  meaningful  stress  and  temperature  is  obtained  by 
taking  either  the  real  or  the  imaginary  part  of  CT,„(x,t)  or  T(x,t),  respectively. 
Substituting  Eqs.  (13)  and  (14)  into  Eq.  (12),  results  in 

-C(to)F  =  r,a(to)o_(i)  (15) 

Since  cy^„(x)  is  given,  in  principle  F(x,cd)  may  be  calculated  using  Eq.  (15). 
For  convenience  V  is  written  as 

F  =  ©  +  z0  (16) 

where  ©  =  ReF  and  0  =  ImF.  Hence,  from  Eq.  (14), 

Re  r  =  2^  +  ©  cos(fflt)  -  0  sin(cot)  (17) 

3.1.1.2  External  Work  Lost 

Using  the  framework  just  developed,  an  expression  is  now  derived  for  the 
external  work  lost  in  terms  of  the  fluctuating  temperature  F.  The  change  in 
temperature  produces  a  thermal  strain.  The  strain  tensor  may  be 
decomposed  into  an  elastic  part  and  a  thermal  part  sj* , 

=  (18) 

Furthermore,  the  thermal  strain  may  be  separated  into  its  spatial  and  temporal 
parts. 


4(x.O  =  £t(x,(o)e'“' 


(19) 
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The  following  identity  derived  by  Brillouin  (1953)  will  be  found  quite  useful 

27t/co  -  - 

^  |[Re(f)Re(/)]rfr=|Re(F/)  =  iRe(F/)  (20) 

where  F  and / represent  two  time-harmonic  complex  qiiantities  with  frequency 
CO,  and  the  overbar  denotes  the  complex  conjugate.  Taking  the  real  values  of 
Eqs.  (13)  and  (19),  substituting  into  Eq.  (2),  and  using  Eq.  (20)  gives  for  Aw, 
the  local  work  lost  per  cycle  per  unit  volmne, 

Aw(x,0)  =  -Jta,„(x)Im(8t(x,co))  (21) 


(The  elastic  part  of  the  strain  field  does  not  contribute  to  the  work  lost  since  it 
is  in  phase  with  the  stress  field).  But 


S/t  (x,  to )  =  a5,„F  =  a6;„  (©  +  i  9) 

(22) 

so  that 

Im(sf,(x,a)))  =  a5,„e 

(23) 

Substituting  Eq.  (23)  into  (21)  gives 

Aw(x,a))  =  -7ia0a„„(x) 

(24) 

Using  Eq.  (3)  the  total  external  work  lost  per  cycle  is  thus 

Ar(©)  =  -7t  J  a  0  (x)  dR 

(25) 

R 


3.1.1.3  Heat  Increment  per  Cycle 

Using  the  framework  developed  for  calculating  the  temperature  field,  an 
expression  is  now  derived  for  the  heat  increment  per  cycle  in  terms  of  the 
fluctuating  temperature  V.  For  an  isotropic  material,  Eq.  (11)  reduces  to 
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(26) 

Taking  the  real  value  of  Eq.  (14),  substituting  into  Eq.  (26)  and  then  using 
Eqs.  (7)  and  (20)  gives  the  local  entropy  produced  per  unit  volume  during  one 
cycle  of  vibration, 

(27) 

Using  Eq.  (8)  the  total  entropy  produced  throughout  the  body  per  cycle  is  thus 

^(®)=-4j*(®..®..+9..e..)‘«  (28) 

^  R 

Since  the  fluctuations  in  temperature  about  7^  are  very  small  (recall  that 
A  TjTg  =  O(l0"^)),  the  heat  increment  per  cycle  LQ  may  be  equated  to  T^AiS. 


3.1.1.4  Equality  of  AW  and  T^AS' 

Using  the  identity 

(29) 

and  the  divergence  theorem,  Eq.  (28)  becomes 


ZAS  = 


n  1 


j  (t©„  @ + w,.  ev.  ds  -  ),.  0 + («., ),.  e)<® 


(30) 


.S  R 

where  S  represents  the  surface  of  the  medium,  and  are  the  components  of 
the  unit  outward  normal  vector  to  the  surface.  Since  the  boundary  of  the 
medium  is  taken  to  be  adiabatic,  F,„,  =  0,„,  +  i 0,„  =  0,  and  therefore 
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Thus,  the  surfacc  term  vanishes  identically,  and  Eq.  (30) 

reduces  to 

J;as = j  ((te,.),.  ® +(te,J.,  9)dR  (31) 

Substituting  Eq.  (16)  into  Eq.  (15)  yields  the  two  real-valued  equations, 

(i®,.L+Cffl6  =  0  (32) 

and 

(^e,J,„-C®©  =  r„CDacT„„(x)  (33) 

Solving  these  two  equations  for  (A:0,„),„  and  {kQ,J,„  and  substituting  into 
Eq.  (31),  yields 

r„A5(a))  =  -7iJaea„„(x)c;i?  (34) 

R 

Comparing  Eq.  (25)  and  Eq.  (34)  gives 

AW:=T„AS  (35) 

Thus,  for  an  isotropic,  homogeneous,  thermoelastic  solid  with  adiabatic 

boundaries,  stressed  time-harmonically,  the  total  work  lost  p>er  cycle  AIT  is 
equal  to  the  heat  increment  per  cycle  T^AS.  This  result  establishes  two 
equivalent  approaches  to  calculating  the  elastothermodynamic  damping:  (1) 

the  mechanical  approach  developed  in  Section  3. 1.1. 2,  and  (2)  the  entropic 
approach  developed  in  Section  3. 1.1. 3. 

3.1.2  A  Counterexample  to  Demonstrate  that  Aw  ^  T^As 

Using  the  framework  developed  in  Section  3.1.1,  a  counterexample  is 
now  constracted  to  demonstrate  that  the  work  lost  per  cycle  per  unit  volume 
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Aw  may  not  be  equal  to  the  local  entropy  produced  per  cycle  per  unit  volxime 
As  times  the  equilibrium  temperature  T^.  Consider  a  homogeneous  isotropic 
Euler-Beroulli  beam  of  rectangular  cross-section  of  thickness  h  with  adiabatic 
boundaries,  subjected  to  an  applied  time-harmonic  bending  moment 
M=M„  ,  as  shown  in  Figure  1 .  The  only  nonzero  stress  is  given  by 

Mz 


=  — 


(36) 


where  /  is  the  second  moment  of  the  cross-sectional  area,  the  z-axis  is  in  the 
thickness  direction,  and  the  x-axis  coincides  with  the  centroidal  axis. 
Substituting  Eq.  (36)  into  the  governing  equation  for  V,  Eq.  (15),  results  in 


C  a  M 

^K-^(iG))F  =  -^^r„(/CD)z 
dz  k  k  I 


—  <z<  — 
2  2 


(37) 


The  adiabatic  boundary  conditions  are 

d 


h 


V  =  0  at  z  =  ±— 
dz  2 


(38) 


Assuming  a  solution  of  the  form 


F  =  ^  ^  sm(n7iz/  h)  (39) 

n=l 

odd 

and  utilizing  the  integral  transform  technique  (Ozisik,  pp.  594-610)  yields 


n-l 


A.= 


aT„M^4h(-Y)^  Q^+iQn^ 
n*+Q,^ 


C  I  Tt^  rP’ 


(40) 


where  a  nondimensional  frequency  has  been  introduced  as  Q  =  (uCh^  j n^k . 
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To  calculate  the  local  work  lost  per  cycle  per  unit  volume,  Eqs.  (39) 
and  (40)  are  substituted  into  Eq.  (24)  yielding 


Aw  = 


(oMMifz  4 


C 


Sc-i) 


odd 


^  ^  ^  (  n\ 

n  +Q  h 


(41) 


To  calculate  the  local  entropy  produced  per  cycle  per  unit  volume,  Eqs.  (39) 
and  (40)  are  substituted  into  Eq.  (27)  resulting  in 


16|.£(-1) 


71 


n=\  m=l 
odd  odd 


nm 


q(q 

(«"*  +  Q^ )( 


^+nV) 


m 


7^ 


cos{mtzf  h)  cos{mTtzf  h) 


(42) 

Motivated  by  the  forms  of  Eqs.  (41)  and  (42),  we  now  introduce  a  normalized 
local  work  lost. 


Aw 

^  —  rp 

and  a  normalized  entropy  produced. 


(43) 


X  = 


ZAs 

(aM,h/iy  ^ 


(44) 


v|/  and  X  plotted  as  functions  of  the  normalized  position,  ^  —  zjh,  in 
Figures  2  and  3,  respectively,  for  tiiree  values  of  the  normalized  frequency  Q. 
For  a  given  frequency,  is  maximum  at  the  outer  surface  of  the  beam, 
^  =  0.5,  whereas  %  is  maximum  at  the  neutral  axis,  ^  =  0.  Thus,  in  general  the 


Normalized  Position,  ^ 


Figure  2.  Normalized  work  lost  \|/  for  an  Euler-Bemoulli  beam  as  a  function 
of  the  normalized  thickness  coordinate  ^  for  three  normalized  frequencies  Q. 


Normalized  Position,  ^ 


Figure  3.  Normalized  entropy  produced  x  for  an  Euler-Bemoulli  beam  as  a 
function  of  the  normalized  thickness  coordinate  ^  for  three  normalized 
frequencies  Q. 
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entropy  is  produced  at  a  different  location  than  where  the  work  is  lost,  i.e. 
Aw 9^7^ As.  As  a  check  on  the  calculations,  the  equality  AW  =7^ AS  was 
verified  explicitly  using  Eqs.  (41)  and  (42). 

To  further  elucidate  the  results  of  this  counterexample,  the  following 
thought  experiment  is  constructed.  Given  a  high-temperature  reservoirs  {?„) 
and  a  low-temperature  reservoir  (7^)  as  shown  in  Figure  4,  a  Camot  heat 
pmnp  supplied  with  an  amount  of  work  AW  =  {l-Tj^fTjj)Q  per  unit  time  is 
used  to  reversibly  transfer  a  quantity  of  heat  {TiJT„)Q  per  unit  time  from  the 
low  temperature  reservoir  and  deliver  a  quantity  of  heat  Q  per  unit  time  to  the 
high  temperature  reservoir  (Callen,  1985).  This  heat  is  then  transferred 
irreversibly  from  the  high  temperature  reservoir  to  the  low  temperature 
reservoir  by  way  of,  say,  a  copper  rod.  Under  steady  state  conditions,  the 
thermodynamic  coordinates  of  both  the  high-temperature  reservoir  and  the 
copper  rod  do  not  change.  However,  an  amount  of  entropy  per  unit  time 
SjN  =  QlTfi  enters  the  copper  rod  from  the  high  temperature  reservoir, 
whereas  an  amovmt  of  entropy  per  unit  time  =  0/7^  leaves  the  copper 
rod  and  enters  the  low  temperature  reservoir.  Hence,  an  amount  of  entropy 
per  unit  time  AS  =  jg/7]_  -  2/7^  is  produced  within  the  copper  rod  and 
transferred  to  the  low-temperature  reservoir  (Alternatively,  this  result  may  be 
readily  obtained  by  using  Eq.  (5)  and  noting  that  the  temperature  varies 
linearly  within  the  copper  rod).  The  result  of  the  entire  process  is  that  the  low 
temperature  reservoir  heats  up  by  an  amount  AQ  =  Tj^AS  =  per 

unit  time,  which  is  precisely  equal  to  AW.  Therefore,  an  amount  of  work  AW 
has  been  totally  dissipated  and  turned  into  heat.  Finally,  note  that  even  though 
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the  work  is  supplied  to  the  Carnot  heat  pump  at  one  spatial  location,  the 
entropy  is  produced  in  the  copper  rod,  an  entirely  different  spatial  location. 

3.1.3  Proof  that  AW  =  T^AS  for  a  Composite  Material 

In  Section  3.1.1  it  was  proven  that  AW=TgAS  for  an  isotropic 
homogeneous  material.  The  equality  AW  =  T^AS  will  now  be  proven  for  a 
composite  material  made  up  of  several  subregions,  each  of  which  may  be 
anisotropic,  and  whose  interfaces  may  be  thermally  imperfect. 

3.1.3.1  Temperature 

Consider  a  composite  with  A^-regions  where  the  yth  subregion  is  labeled 
and  has  siuface  S-'  as  shown  in  Figure  5.  The  thermal  conductivity  and  the 
linear  thermal  expansion  are  now  described  by  the  tensors,  and  a/„, 
respectively  (There  is  implied  summation  only  on  repeated  subscripts).  The 
temperature  field  in  each  subregion  is  governed  by  Eq.  (10).  Let  «/  be  the 
unit  outward  normal  vector  to  subregion  R\  The  boundary  and  interface 
conditions  are  as  follows  (Ozisik,  pp.  594-596).  The  outer  boundary  of  the 
composite  is  taken  to  be  adiabatic: 

(9,n,y=0  (45) 

(only  for  that  part  of  the  yth  subregion  that  constitutes  a  part  of  the  outer 
boundary  of  the  composite).  The  interface  conditions  are:  (1)  the  continuity 
of  the  heat  flux  vector 

{q,n, )'  +  {q,n,y  =  0 ,  «[  =  -nj  j^r 


(46) 
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and  (2)  the  condition  of  imperfect  thermal  contact 

=f‘’(T‘-r),  J*r  (47) 

(only  for  the  common  boimdary  of  the  subregions  and  R’'),  where 
is  the  thermal  contact  conductance  between  adjacent  subregions  R^ 
and  R'.  Eq.  (47)  represents  the  discontinuity  of  temperature  between 
subregions  i?-'  and  R"".  For  a  perfect  thermal  interface  oo,  and  Eq.  (47) 
reduces  to  the  familiar  condition  of  continuity  of  temperature,  V  =  f. 

As  an  aside  note  that  the  concept  of  a  thermal  contact  conductance  may 
be  used  not  only  to  model  a  distinct  imperfect  thermal  interface  but  also  an 
interphase  layer  which  is  sufficiently  thin  so  that  a  linear  temperamre 
distribution  may  be  assumed.  Then / may  be  thought  of  as  the  quotient  of  two 
parameters:  (1)  an  effective  interphase  thickness  h,  and  (2)  an  effective 
interphase  thermal  conductivity  k,  such  that  /  =  k//z.  Combining  the 
interface  conditions,  Eq.  (46)  and  Eq.  (47),  gives 

qinj=-K^’'^  J  =q^ni  (48) 

From  this  equation  it  can  be  clearly  seen  that  the  continuity  of  heat  is 
maintained  across  an  interphase  layer  whereas  the  temperature  is 
discontinuous  due  to  the  finite  distance  and  finite  thermal  conductivity  of  the 
inteiphase  layer. 

Substituting  Eqs.  (13)  and  (14)  into  Eq.  (10)  results  in 


(49) 
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3.L3.2  External  Work  Lost 

The  external  work  lost  is  calculated  using  the  steps  presented  in 
Section  3.1. 1.2.  Eqs.  (18),  (19),  and  (21)  remain  the  same  except  now  each  of 
the  quantities  is  decorated  with  a  superscript  (  y  for  the  subregion  R\  Since 
the  coefficient  of  thermal  expansion  is  now  a  tensor,  in  place  of  Eq.  (22) 

Sto*(x,o>)  =  a/„F  =a/„(©-'+ieO  (50) 

so  that 

lm(s/^*(x,0)))  =  ai,e^  (51) 

In  place  of  Eq.  (24) 

Aw^(x,co)  =  -7ta/„(x)ai,e^  (52) 

Finally,  the  total  work  lost  per  cycle  in  the  composite  is 

AIF(o))  =  -7r2]Jai,(x)a/„0^<//?  (53) 

3.1.3.3  Heat  Increment  per  Cycle 

Taking  the  real  value  of  Eq.  (14),  substituting  into  Eq.  (1 1)  and  then  using 
Eqs.  (7)  and  (20)  gives 

(54) 

The  entropy  produced  per  unit  time  at  the  thermally  imperfect  interfaces  due 
to  a  discontinuity  in  temperature  is  calculated  as  follows.  The  change  in 
entropy  per  unit  time  of  a  subregion  by  heat  flux  through  a  surface  is 
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-\  {qjnj IP)dS.  Given  any  two  subregions  and  R' ,  the  entropy 
produced  per  unit  time,  A 5^'',  at  their  mutual  interface,  S^'  =  is  thus 


Therefore,  the  entropy  produced  per  unit  time  at  all  the  interfaces  is  given  by 
the  following  summation  over  the  subregions 

^^interface  =  J ^dS  (56) 

y=i  ij  T-' 

Consistent  with  the  approximations  of  the  linear  theory  of 
elastothermodynamic  damping,  the  factor  l/  P  is  approximated  with' 


X  = Y  =  — fl-Re(FV“')/2;l  (57) 

r  r„+Re(FV“')  7;L  ^  "J 

Substituting  Eqs.  (6)  and  (57)  into  Eq.  (56),  using  Eq.  (14)  and  integrating 
over  one  cycle  with  the  aid  of  Eq.  (20)  gives  the  entropy  produced  per  cycle  at 
the  thermal  interfaces  due  to  a  discontinuity  in  temperature 


■7T  "  f 

^interface  f 

The  first  summation  may  be  rewritten  as  a  sum  over  the  interfaces.  Using  the 
continuity  of  heat  flux  condition  Eq.  (46)  and  the  fact  that  «/  =  for 
adjacent  subregions,  the  first  summation  is  identically  zero,  hence 


ULni{v.i‘"‘-ny'/Tj\dS  (58) 
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Substituting  Eqs.  (54)  and  (59)  into  Eq.  (8)  gives  the  total  entropy  produced 
per  cycle  throughout  the  composite 


(60) 


As  in  the  isotropic  case  the  heat  increment  per  cycle  IsQ  may  be  equated  to 
TAS. 

3.1.3.4  Equality  of  AfV  and  T^AS 
Using  the  identity 

•  ®  ©)./  .m  )>/  ®  (61) 

and  the  divergence  theorem  in  Eq.  (60),  the  surface  terms  vanish  identically 
yielding 

+fee4)„e^]<«  (62) 

Substituting  Eq.  (16)  into  Eq.  (49)  yields  the  real- valued  equations, 

{kL&!.)„*CW  =0  (63) 

and 

te9;J.,-C>CDe^  =  r,coai,o/.(i)  (64) 

Solving  Eq.  (63)  for  (ki,@Q,i  and  Eq.  (64)  for  (*,1,64).;  and  substituting  into 
Eq.  (62)  yields 

T„AS{p)  =  -TT^  j  a/„(x)ai,e"  dR 


(65) 
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Comparing  Eq.  (53)  to  Eq.  (65)  gives 

AW=T,AS  (66) 

Thus,  for  a  thermoelastic  composite  material  consisting  of  anisotropic 
subregions  with  imperfect  thermal  interfaces,  with  an  adiabatic  outer 
boundary,  and  loaded  time-harmonically,  the  total  work  lost  per  cycle  AfF  is 
equal  to  the  total  heat  increment  per  cycle  7^  AS .  This  result  establishes  two 
approaches  to  calculating  the  elastothermodynamic  damping  in  composite 
materials:  (1)  the  mechanical  approach  developed  in  Section  3. 1.3. 2,  and  (2) 
the  entropic  approach  developed  in  Section  3. 1.3. 3. 


3.2  SOLUTION  OF  THE  INHOMOGENEOUS  HEAT-CONDUCTION 
PROBLEM  FOR  A  GENERAL  COMPOSITE  MEDIUM  USING  THE 
INTEGRAL-TRANSFORM  TECHNIQUE 

In  order  to  use  either  the  mechanical  approach  or  the  entropic  approach  to 
calculate  the  elastothermodynamic  damping  in  a  composite  material,  the 
temperature  field,  governed  by  Eq.  (10),  must  first  be  determined.  The  first 
step  in  this  solution  is  to  view  the  inhomogeneous  term  in  Eq.  (10)  due  to  the 
thermoelastic  effect  as  a  heat-generation  term.  Ozisik  (1980,  pp.  594-610)  has 
solved  the  inhomogeneous  heat  conduction  problem  (the  heat  conduction 
problem  including  heat  generation)  for  a  composite  medium  consisting  of 
isotropic  subregions  with  thermally  imperfect  interfaces  using  the  integral 
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transform  technique.  In  this  section  the  solution  given  by  Ozisik  is  extended 
to  a  composite  medium  consisting  of  anisotropic  subregions. 

The  basic  steps  in  the  solution  of  the  inhomogeneous  heat-conduction 
problem  using  the  integral-transform  technique  are  as  follows. 

(1)  The  integral  transform  and  inversion  formulae  are  developed  using  the 
eigenfunctions  of  the  associated  homogeneous  problem. 

(2)  By  application  of  the  integral  transform,  the  partial  differential  heat- 
conduction  equation  in  space  and  time  is  reduced  to  an  ordinary  differential 
equation  in  time. 

(3)  The  resulting  ordinary  differential  equation  is  solved  subject  to  the 
transformed  initial  condition. 

(4)  The  transform  of  the  temperature  is  inverted  using  the  inversion  formula 
to  obtain  the  desired  solution. 

3.2.1  The  Time-Dependent  Inhomogeneous  Heat-Conduction  Problem 
for  a  Composite  Medium 

Consider  a  composite  medium  consisting  of  iV-anisotropic  subregions 
with  heat  generation.  As  in  Section  3.1.3  the  jUi  subregion  is  labeled  R\  has 
surface  5'-',  and  a  unit  outward  normal  vector  «/  as  shown  in  Figure  5.  The 
temperature  field  in  each  subregion  is  governed  by 

^ter.i).,+^r'(*,0=|:r'(x.o  r>o  (67) 
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where  g^  (x,t)  is  the  heat  generation  per  unit  volume  in  subregion  The 
boundary  and  interface  conditions  are  identical  to  those  described  in 
Section  3 . 1 .3 . 1 .  The  initial  condition  is  taken  to  be 

P(x,t)  =  T„  at  t  =  0  (68) 

3.2.2  Development  of  the  Integral  Transform  and  Inversion  Formulae 

The  first  step  in  solving  this  inhomogeneous  time-dependent  heat- 
conduction  problem  by  the  method  of  integral  transform  is  to  develop  the 
integral  transform  and  inversion  formulae.  At  the  foundation  of  this 
development  is  the  solution  of  the  homogeneous  problem.  In  order  to  obtain 
the  homogeneous  solution,  the  classical  technique  of  separation  of  variables  is 
used. 

3.2.2.1  Solution  of  the  Homogeneous  Heat-Conduction  Problem 

Assuming  a  temperature  solution  of  the  form  P(x,t)  =  ^-'(x)r(t),  the 
homogeneous  version  of  Eq.  (67)  may  be  separated  into  spatial  and  temporal 
equations, 

r(0  =  e-^''  (69) 

and 

PO) 

where  P  is  the  separation  constant.  The  boundary  and  interface  conditions 
become 


(71) 
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for  the  outer  boundary,  and 

-w.«/=/"(v-f)=-o.;x'  (72) 

for  the  common  boundary  of  the  subregions  and  R' .  Eqs.  (70),  (71), 
and  (72)  constitute  an  eigenvalue  problem  for  the  eigenfunctions  «t»-'(P„,x)  in 
terms  of  the  eigenvalues  P„  (Unless  otherwise  stated,  only  the  subscripts  /  and 
m  are  tensorial  indices.). 

3.2.2.2  Derivation  of  the  Orthogonality  Condition 

The  eigenfunctions  obtained  by  solving  the  eigenvalue  problem  have  a 
very  special  property,  a  generalized  orthogonality.  This  generalized 
orthogonality  property  is  the  heart  of  the  integral  transform  technique.  In 
order  to  derive  the  orthogonality  condition,  Eq.  (70)  is  written  for  two 
different  eigenfunctions  <|)''(P„,x)  s  and  (|)''(P^,x)  s  as 

(73) 

and 

(74) 

Eq.  (73)  is  multiplied  by  and  Eq.  (74)  by  the  resulting  equations  are 
subtracted  and  integrated  over  the  subregion  R^ .  Summing  over  all 
subregions  gives 

(p"  J r'^'dR  =  i  ¥’]dR  (75) 

7=1  Ri  7=1  Ri 

Using  the  identity  of  Eq.  (61)  and  the  divergence  theorem,  Eq.  (75)  becomes 
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(p; V'  - to,i'y'']»/  (76) 

y=i  M  sj 

where  the  volume  integrals  have  canceled.  Using  the  outer  boundary  and 
interface  conditions,  Eqs.  (71)  and  (72),  and  the  fact  that  nj  =  -n'  for  adjacent 
subregions  and  J?'',  the  surface  integrals  in  Eq.  (76)  vanish  identically. 
Thus, 


7=1  RJ 

which  yields  the  orthogonality  relation 


>1 


n^p 
n  =  p 


where 


7=1  RJ 


(77) 


(78) 


(79) 


3^.2.3  The  Integral  Transform  and  Inversion  Formulae 

Now  that  the  orthogonality  relation  for  the  eigenfunctions  has  been 
obtained,  the  integral  transform  and  inversion  formulae  may  be  derived. 
Consider  the  following  representation  of  a  function  P(x,t),  defined  in  a 
finite  region  i?-',  in  terms  of  the  eigenfunctions  (|)^(p„,x)  of  the  previous 
eigenvalue  problem, 

r(x,/)=i;4.(o<i''(p..*) 

IJ=1 


(80) 
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where  the  summation  is  over  all  eigenvalues  Multiplying  Eq.  (80)  by 
integrating  over  region  R\  and  summing  over  all  subregions  gives 

=  (81) 

>1  kJ  n=l  j=\  gJ 

Using  the  orthogonality  relation,  Eq.  (78),  and  solving  for  A,,  (t)  yields 


w = Z  I  f 

J=l 


(82) 


Eq.  (82)  is  substituted  into  Eq.  (80)  and  the  resulting  expression  is  split  into 
two  parts  to  define  the  desired  integral  transform  and  inversion  formulae 

^  .  r  ■ 

Integral  transform:  f  (P„ , t)  =  ^  j  <|)-'  (p„  ,x)  V  (x, t) dR  (83) 

>1  Ri 

Inversion  formula:  P(\,t)  =  ^(P«  ’0  (84) 


3.2.3  Integral  Transform  of  the  Inhomogeneous  Heat-Conduction 
Problem 

Having  established  the  integral  transform  and  inversion  formulae,  the 
inhomogeneous  conduction-problem  stated  in  Section  3.2.1  may  now  be 
solved.  The  spatial  derivatives  in  Eq.  (67)  may  be  removed  with  the 
application  of  the  integral  transform.  Multiplying  Eq.  (67)  by  C^cj)-'’", 
integrating  over  subregion  R\  and  slimming  over  all  subregions  gives 


2  fV''(iir.:).,rfs+|;  'g\x,/)rfR =42'^' 

7=1  Rj  7=1  Ri  7=1  Ri 
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Using  the  integral  transform  formula  Eq.  (83),  the  summation  on  the  right  is 
simply  f{t) .  Using  the  identity  given  in  Eq.  (61)  and  the  divergence  theorem 
twice,  the  first  summation  in  Eq.  (85)  becomes 

i;  =  I; \l>HKTinl]dS-f^\Tikll,r  ni)dS 

j=\Rj  j=isJ  ^  J=^sJ  (35) 

*'ZiTik/jr),.dR 

J=^  Rj 

Using  the  outer  boundary  and  interface  conditions,  Eqs.  (71)  and  (72),  and  the 
fact  that  «/  =  -Hi  for  adjacent  subregions  and  R'',  the  surface  integrals  in 
Eq.  (86)  vanish  identically.  Also,  from  Eq.  (70),  =  ,  so 

that  Eq.  (86)  reduces  to 

(87) 

where  the  integral  transform  formula  Eq.  (83)  has  been  used.  Substituting 
Eq.  (87)  into  Eq.  (85)  yields 

^ni)+Km=m  (ss) 

at 

where  an  integral  transform  of  g^{x,t)  has  been  defined  as 

=  (89) 

>1 

3.2.4  Solution  for  the  Integral  Transform  of  the  Temperature  and 
Inversion 

The  integral  transform  of  the  initial  condition  Eq.  (68)  is 
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-/=1  Rj 

The  solution  of  Eq.  (88)  subject  to  this  initial  condition  is 


(90) 


T{t)  =  e 


- 


To  + 


Je'’""^(T)i/T 


(91) 


Using  the  inversion  formula  Eq.  (84),  the  temperature  field  of  the  composite 
medium  is  given  by 


f„+^e^"^g(x)dx 


(92) 


Since  this  solution  incorporates  adiabatic  outer  boundaries,  the  first 
eigenvalue  is  Pi=0  and  the  first  eigenfunction  is  With  this 

information  it  may  be  shown  that 

e-*'-' =  (93) 

Hence,  Eq.  (92)  reduces  to 


=  r.  A  (94) 

w=l  n  o 

This  equation  represents  the  temperature  field  of  a  composite  medium 
consisting  of  anisotropic  subregions  with  arbitrary  heat  generation  and  whose 
interfaces  may  be  thermally  imperfect. 
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3.3  ELASTOTHERMODYNAMIC  DAMPING  IN  A  GENERAL 
COMPOSITE  MEDIUM 

Now  that  an  expression  for  the  temperature  field  has  been  obtained  for  a 
composite  material  with  arbitrary  heat  generation,  expressions  may  now  be 
derived  for  the  temperature  field,  local  work  lost,  local  entropy  produced,  and 
the  elastothermodynamic  damping.  Comparing  Eq.  (10)  with  Eq.  (67),  it  is 
clear  that  for  the  problem  of  elastothermodynamic  damping 

g^(x,0  =  -W„^  (95) 

ot 

Substituting  Eq.  (13)  into  Eq.  (95)  gives 

g^(x,0  =  -TSm)ai  (Ti,(x)c'“'  (96) 

Substituting  Eq.  (96)  into  Eq.  (89)  results  in 

m  =  (97) 

where 

(98) 

Substituting  Eq.  (97)  into  Eq.  (94)  and  carrying  out  the  integration  in  time 
yields  the  following  expression  for  the  temperature  field  of  the  composite 

P(x,t)  =  -«■'■')  (99) 

The  first  term  in  parenthesis  represents  the  time-harmonic  solution  whereas 
the  second  term  represents  the  transient  solution.  Herein,  only  the  time- 
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harmonic  solution  is  of  interest.  Comparing  Eq.  (99)  with  Eq.  (14),  it  is  clear 
that  the  fluctuating  temperature  is  given  by 


+  N„ 


(100) 


Substituting  Eq.  (100)  into  Eq.  (52)  gives  the  local  work  lost  per  unit 


volume  per  cycle, 


A«^(x,c»)  =  „r,|;-^ioi,oi,(x).l."(x)  (101) 

+P„  K 

Substituting  Eq.  (100)  into  Eq.  (54)  gives  the  local  entropy  produced  per  unit 
volume  per  cycle, 


”  "  cofco^+B^P^)  r  L  .  . 

(102) 


'  '  iv,  . . 

As  shown  in  Section  3.1.3,  the  total  external  work  lost  over  the  entire 


composite  per  cycle  may  be  calculated  using  either  the  local  work  lost  or  the 
local  entropy  produced.  From  the  form  of  Eqs.  (101)  and  (102)  it  is 
mathematically  much  easier  to  use  the  local  work  lost  expressions.  Thus, 
substituting  Eq.  (100)  into  Eq.  (53),  the  total  external  work  lost  per  cycle  over 
the  entire  composite  is  given  by 


AFF(co)  =  7i2;X 

n=\ 


«>P^  4 


(103) 


The  elastothermodynamic  damping  of  a  composite  medium  consisting  of  an 
arbitrary  number  of  subregions,  each  of  which  may  be  anisotropic  and  whose 
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thermal  interfaces  may  be  imperfect,  subjected  to  any  time-harmonic  stress 
field,  may  be  calculated  using  Eq.  (103). 


3.4  ELASTOTHERMODYNAMIC  DAMPING  IN  LAYERED 
COMPOSITES 

The  results  of  Section  3.3  are  now  specialized  to  a  composite  material 
consisting  of  N  isotropic  layers  in  a  rectangular,  cylindrical,  and  spherical 
coordinate  system  (i.e.  an  iV-layer  slab,  cylinder,  and  sphere,  respectively), 
subjected  to  any  stress  field  so  long  as  the  resulting  heat  conduction  can  be 
described  by  a  single  spatial  coordinate  orthogonal  to  the  layering  (see 
Figure  6).  Many  problems  of  current  engineering  interest  can  be  modeled  as 
special  cases  of  the  JV-layer  solution.  (1)  Cartesian  coordinate  system: 
laminated  composites  in  bending  and  extension;  aligned  platelet  composite;  or 
a  periodic  array  of  labs.  (2)  Cylindrical  coordinate  system:  a  unidirectional 
fiber-reinforced  composite  where  each  fiber  may  consist  of  several  concentric 
cylinders.  (3)  Spherical  coordinated  system:  a  particulate  composite  where 
each  particle  may  be  solid  or  hollow,  or  consist  of  several  concentric  spheres, 
for  example,  an  aluminum  matrix  filled  with  hollow  AI2O3  microspheres. 
Specific  results  are  given  for  the  following  problems.  (I)  A  periodic  array  of 
slabs  (rectangular  symmetry)  with  a  two-layer  representative  volume  element 
subjected  to  the  following  mechanical  stress  states:  (a)  a  time-harmonic 
uniform  stress  perpendicular  to  the  layering,  and  (b)  a  time-harmonic  uniform 
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strain  parallel  to  the  layering.  (II)  A  symmetric  three-layer  plate  in  pure 
bending  with  a  stress  field  obtained  from  the  classical  laminate  theory.  (Ill)  A 
hollow  two-layer  concentric  sphere  subjected  to  a  uniform  time-harmonic 
radial  load. 

The  fluctuating  temperature  field  of  a  general  composite  medium 
subjected  to  any  time-harmonic  loading  is  given  by  Eq.  (100).  For  an  iV-layer 
composite  with  isotropic  layers,  Eq.  ( 1 00)  reduces  to 


1(0 


(104) 


where  x  is  the  spatial  coordinate  normal  to  the  layering  as  shown  in  Figure  6. 
(In  this  section  a  particular  layer  will  be  referenced  with  a  subscript  instead  of 
a  superscript  as  in  Section  3.3.)  But  now  Eq.  (79)  reduces  to 


=  (105) 

where  is  the  Sturm-Liouville  weight  function  (p  =  0,1,2  for  rectangular, 
cylindrical,  and  spherical  symmetry,  respectively).  Also,  Eq.  (98)  reduces  to 

(106) 

J-^  X . 

The  local  work  lost  per  unit  volume  per  cycle  in  a  general  composite 
medium  is  described  by  Eq.  (101).  For  an  iV-layer  composite  with  isotropic 
layers,  this  expression  reduces  to 


Awj  (x,©)  =  nT„a.jO^  (x)^] 


top'  4 


(107) 
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The  local  entropy  produced  per  unit  volume  per  cycle  in  a  general  composite 
is  described  by  Eq.  (102).  For  an  iV-layer  composite  this  equation  reduces  to 


A.  (X  y 


(108) 


The  total  external  work  lost  per  cycle  for  a  general  composite  medium  is 
described  by  Eq.  (103).  This  expression  has  the  same  form  for  an  Allayer 
composite  but  now  N„  and  L„  are  given  by  Eqs.  (105)  and  (106),  respectively. 


3.4.1  Elastothermodynamic  Damping  in  a  Periodic  Array  of  Slabs 

Consider  a  periodic  array  of  slabs  of  two  isotropic  materials  with 
imperfect  thermal  interfaces  as  shown  in  Figure  7.  A  unit  cell  is  shown  in  an 
inset  in  the  figure  and  consists  of  a  layer  of  material  (1)  with  thickness  a  and  a 
layer  of  material  (2)  with  thickness  (b-a).  Unless  otherwise  stated,  in  the 
following  the  subscripts  (  ),  and  (  )2  will  refer  to  materials  (1)  and  (2), 
respectively.  The  planes  x  =  0  and  x  =  b  are  adiabatic  due  to  reflection 
symmetry  on  these  planes.  The  eigenfunctions  for  each  layer  can  be 
determined  using  the  method  of  Ozisik  (1980,  pp.  294-334),  giving 

(|),„(x)  =  cos(Y„x/a) 

K  (x)  =  A  cos(ri„x/b)  -i-  B„  sin(Ti„x/b) 


(109) 

(110) 
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where 


4,= 

B  = 


k 

cos(Y„)--^Y„sin(y„) 

k 

cos(Y„)--ry„sin(Y„) 


cos(Ti„a/6)  +  K  sin(Y„ )  sin(Ti„a/h) 
sin(Ti„a/h)-  ^sin(Y„)  cos(Ti„a/6) 


(111) 


and 


b-a' 


(112) 


z,=(C/k\a\  z,=(C/k\ib-af .  K  =  ^(kC}J{kC\ 

The  constant  x  is  a  characteristic  time  of  heat  conduction,  and /is  the  thermal 
contact  conductance  for  the  interface.  The  transcendental  equation  for  the 
eigenvalues  is 


k 

-7-Y«sin(Y„)-cos(Y„)  cos(Ti„a/6)  sm{\\„alb) 
fa 

/:sin(Y„)  -sin(Ti„a/6)  cos(t|„^7/Z?) 

0  -sin(qj  cos(ti„) 


=  0 


(113) 


The  following  two  canonical  time-harmonic  stress  states  are  considered: 
(1)  a  uniform  uniaxial  stress  perpendicular  to  the  layers 


^xtl  ~  ^xit2  ~  ^ 


otherwise 


(114) 


(2)  a  uniform  strain  parallel  to  the  layers 


^yyi  ® 
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so  that 


cyyy,  =  £,8  and 


^Ui  ~  ^iJz  ~  ^  Otherwise 


(115) 


where  cy  and  s  are  constants,  E  is  the  Young's  modulus,  and  the  time-harmonic 
factor  e""'  is  implied.  These  two  stress  states  will  be  referred  to  as 
problem  (1)  and  problem  (2),  respectively. 

Using  Eq.  (104)  the  fluctuating  temperature  in  layers  1  and  2  for  both 
problems  (1)  and  (2)  is  given  by 


=  ;  =  1,2  (116) 
»=1  +yn 

where  a  dimensionless  frequency  Q,  =  (Dt,  has  been  defined,  and 


n  _  6  '  a.  ' 

— - 7;; - 

-I  +^/ 

h‘^  c,' 

where  F  =  1  for  problem  (1)  and  F  =  E^jE^  for  problem  (2),  and 

/,=sin(Y„)/Y„ 

h  =— {4[sin(Ti„)-sin(q„a/6)]-F„[cos(Ti„)-cos(Ti„a/i)]} 

^  =  ^+7^sin(2Y„) 

2  4y„ 


(117) 
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h  =|(^l-04'+4')+;^(4-4')[sin(2Ti„)-sin(2Ti„a/Z))] 

-  A^n  [cos(2ti„  )  -  cos(2Ti„a/6)] 

2ti„ 

(118) 

Now  intuitively  as  Q,  0  the  time  period  goes  to  infinity.  During  each 
cycle  heat  has  plenty  of  time  to  conduct  between  'cold'  and  'hot'  slabs,  and 
essentially  isothermal  conditions  are  obtained.  In  the  other  extreme  as 
Q,  ->  00  the  time  period  goes  to  zero,  and  heat  has  very  little  time  to  conduct. 
Therefore,  with  the  exception  of  the  interface  where  there  will  always  be  some 
heat  conduction  no  matter  how  high  the  frequency,  essentially  adiabatic 
conditions  are  obtained.  These  intuitive  observations  will  now  be  used  as  a 
check  on  Eq.  (116). 

Suppose  the  unit  cell  is  initially  at  rest  with  equilibrium  temperature  7^ . 
At  time  /  =  0  the  stresses  (Tjyt,  and  are  suddenly  applied  to  layers  (1) 
and  (2),  respectively,  and  held  constant  thereafter.  Within  the  linear  theory  of 
elastotheraiodynamic  damping,  fi’om  Eq.  (1)  the  instantaneous  changes  in 
temperature  of  the  two  layers  due  to  the  thermoelastic  effect  are 

ATj=-T,{ac^/Cl  7  =  1,2  (119) 

Upon  the  completion  of  the  heat  conduction  process,  the  unit  cell  will  be  at  a 
homogeneous  equilibrium  temperature  given  by 

{a/b){a(y^ ),+(!-  alb){a<y^  \ 

J  0  °  (a/h)C,+(l-a/Z>)C2 


7  =  1,2  (120) 
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It  has  been  verified  that  Eq.  (116)  correctly  reduces  to  this  isothermal  limit  as 
Qj-»0.  Under  strictly  adiabatic  conditions  (i.e.  when  the  thermal 
conductivity  k  =  0),  the  temperature  in  each  layer  is  given  merely  by  the 
thermoelastic  effect,  Eq.  (119);  that  is 

V,=-Tiax!jC)^  j  =  \,2  (121) 

Once  again,  with  the  exception  of  the  immediate  vicinity  of  the  interface,  it 
was  shown  that  Eq.  (116)  correctly  reduces  to  this  adiabatic  limit  as  O,  — >  oo. 
In  passing  note  that  the  parameters  a/C  for  problem  (1)  and  EaJC  for 
problem  (2)  are  the  most  significant  factors  determining  the  magnitude  of  V. 
The  numerical  values  of  these  parameters  for  some  pertinent  engineering 
materials  are  presented  in  Table  1. 


Table  1.  Thermoelastic  constants  and  dimensionless  parameters 


E 

(GPa) 

k 

(J/smK) 

a 

(10-6) 

(1/K) 

C 

(106) 

(J/m^K) 

a/C 

(10-12) 

(J/m3) 

£a/C 

^0 

A1 

70 

111 

24 

2.43 

9.71 

0.68 

0.0302 

Steel 

200 

52 

12 

3.82 

3.14 

0.63 

0.0147 

Zn 

103 

113 

40 

2.73 

14.5 

1.5 

0.1121 

Mg 

44 

154 

27 

1.79 

15.1 

0.67 

0.0340 

Ti 

116 

39 

8.4 

2.34 

3.59 

0.42 

0.0066 

AI2O3 

350 

29 

9 

3.19 

2.82 

0.99 

0.0167 

SiC 

460 

90 

4.3 

4.34 

0.99 

0.46 

0.0037 
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Numerical  results  for  problem  (1)  for  a  SiC-Al  periodic  layered 
composite  will  now  be  presented.  Let  layer  (1)  be  SiC  and  layer  (2)  be  Al. 
The  volume  fraction  of  the  SiC  "inclusions"  is  denoted  by  =  ajb .  A 
normalized  fluctuating  temperature  ©  is  introduced  as 


0  = 


(122) 


(T.'y.alcX 

In  polar  form  0  =  ||0|[  exp(/  O) ,  where  ||0||  is  the  magnitude  and  O  is  the  phase 
of  0.  Consider  first  a  thermally  perfect  interface.  In  Figures  8  and  9  ||0||  and 
O  are  plotted,  respectively,  verses  the  normalized  position  ^  =  xjb  with 
=  ^sic  ^  ^  parameter.  Note  that  for  Qgjc  « 1  the  conditions  are  essentially 
isothermal  and  Osti.  As  Qgic  increases  significant  temperature  gradients 
develop  and  become  localized  near  the  interface.  At  very  high  frequencies, 
Qsic  the  conditions  are  essentially  adiabatic  and,  once  again,  Ostc 

(with  the  obvious  exception  of  the  immediate  vicinity  of  the  interface  as 
discussed  previously). 

Now  consider  a  thermally  imperfect  interface.  A  normalized  thermal 
contact  conductance  ^  is  introduced  as 


c 


K 


(123) 
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adiabatic  plane  Interfece  adiabatic  plane 


Noraialized  Position  ^ 

Figure  8.  Magnitude  of  the  normalized  fluctuating  temperature,  ||0||,  as  a 
function  of  the  normalized  position  ^  for  problem  (1)  with  a  perfect  thermal 
interface,  for  the  composite  SiC-Al,  for  foxor  normalized  frequencies  QgiC’ 
for  a  SiC  volume  flection  of  )^  =  0.5. 
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adiabatic  plane  ,  j-  u  i 

^  mtertace  adiabatic  plane 


Noraialized  Position,  ^ 


Figure  9.  Phase  of  the  fluctuating  temperature,  O,  as  a  function  of  the 
normalized  position  §  for  problem  (1)  with  a  perfect  thermal  interface,  for  the 
composite  SiC-Al,  for  four  normalized  frequencies  QjiC’  ^  volume 
fraction  of  =  0.5. 
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Moreover,  if  the  thermal  contact  conductance  is  used  to  model  the  presence  of 
an  extremely  thin  interphase  layer  of  thickness  h  and  thermal  conductivity  k 
(so  that  a  linear  temperature  variation  across  the  thickness  may  be  assumed), 
then  /  =  k//2  and  therefore 


K 


K  a 

T,h 


(124) 


In  the  limit  as  /z  0  or  k  oo,  /  oo  and  the  conditions  corresponding  to  a 
thermally  perfect  interface  are  recovered. 

In  Figure  10  ||0||  is  plotted  versus  the  normalized  spatial  coordinate  ^  with 
^  as  a  parameter  for  a  normalized  j&equency  of  Qgic  Ot  will  be  shown 
later  that  this  particular  frequency  corresponds  to  a  maximimi  in  the  total 
damping.)  A  comparison  of  Figures  8  and  10  reveals  that  conditions 
corresponding  to  a  thermally  perfect  interface  are  obtained  as  ^  approaches 
100.  Conversely,  as  ^  approaches  0.1  the  two  layers  become  thermally 
isolated.  For  a  typical  interphase  layer  k  is  of  the  order  of  k^  and 
h/a  =  o{l0~^).  Therefore,  from  Eq.  (124),  ^  is  of  the  order  of  100,  i.e.  a 
condition  corresponding  to  a  thermally  perfect  interface.  One  may  conclude 
that  a  typical  interphase  layer  will  have  little  effect  on  the  temperature  and 
thus  on  the  elastothermodynamic  damping.  This  fact  will  be  further 
elucidated  in  the  following  analysis. 

Using  Eq.  (107)  the  local  work  lost  per  cycle  in  each  layer  Aw^(x,(»), 
j  =  1,2,  is  given  by 


(x,  to)  =  7rr„  (a^^a) ,  (acr**  /C),  ^ 

W=1 


(125) 
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Normalized  Position  ^ 


Figure  10.  Magnitude  of  the  normalized  fluctuating  temperature,  ||0||,  as  a 
function  of  the  normalized  position  4  for  problem  (1)  with  an  imperfect 
thermal  interface,  for  the  composite  SiC-Al,  for  fom  normalized  thermal 
contact  conductances  for  a  normalized  frequency  Qgic  =  10,  and  for  a  SiC 
volume  fraction  of  Vj  =  0.5. 
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A  normalized  local  work  lost  \\f  is  introduced  as 


A  physical  feel  for  the  denominator  may  be  developed  as  follows.  A  rod  of 
one  unit  volume  at  an  equilibrium  temperature  is  subjected  to  a  uniaxial 
tension  of  magnitude  t  under  adiabatic  conditions;  then  V  =  -xTgafC.  While 
holding  T  fixed,  the  rod  is  allowed  to  absorb  heat  from  a  reservoir  at  until  it 
returns  to  its  equilibrium  temperature  7^.  Next,  the  rod  is  unloaded 
adiabatically;  then  V  =  +xTj(xlC.  Once  again,  the  rod  is  allowed  to  give  up 
heat  to  a  reservoir  at  until  the  rod  temperature  becomes  2^ .  The  total  work 
lost  during  this  cycle  is  JC.  By  way  of  example,  for  steel, 

a  =  12x10"* and  C  =  3.8x10*  If  we  choose  t  =  100 MPa, 

Tg  =  300.K^,  and  k  =  1  (plane  stress),  then  Aw  =  114  .  Interestingly,  the 

temperature  increase  of  the  rod  due  to  this  conversion  of  work  into  internal 
heat  energy  is  only  AT  =3x10“* AT.  For  use  in  the  next  section,  the  total 
entropy  produced  during  this  cycle  is  x^a^jC.  For  the  example  under 
consideration.  A?  =  0. 3  8  Jlm^K . 

For  the  case  of  a  thermally  perfect  interface,  vj/  is  plotted  in  Figure  1 1 
versus  the  normalized  position  ^  with  the  normalized  frequency  flgc  as  a 
parameter.  The  maximum  local  work  lost  is  in  the  aluminum  layer  for  all  the 
frequencies  shown.  Note  that  the  local  work  lost  is  actually  negative  in  the 
SiC  layer.  At  first  sight  this  may  appear  disturbing  for  it  implies  that  work  is 
being  produced  rather  than  converted  into  heat.  However,  the  only  constraint 
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Normalized  Position,  4 


Figure  11.  Normalized  work  lost  \|/  as  a  function  of  the  normalized  position  ^ 
for  problem  (1)  with  a  perfect  thermal  interface,  for  the  composite  SiC-Al,  for 
four  values  of  the  normalized  frequency  QsiC’  ^  volume  fraction 

ofF^  =0.5. 
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imposed  by  the  second  law  of  thermodynamics  is  that  the  total  entropy 
produced  throughout  the  composite  (a  closed  thermodynamic  system  which 
does  not  exchange  heat  with  its  surroundings)  during  one  cycle  of  loading  is 
either  zero  or  positive,  i.e.  A<S'>0.  It  follows  from  Eq.  (66)  that 
=  AW  >  0,  i.e.,  the  total  work  lost  throughout  the  composite  during  one 
cycle  should  be  zero  or  positive.  It  follows  that  'P  >  0.  From  an  examination 
of  Figure  15,  it  is  apparent  that  this  second-law  constraint  is  clearly  satisfied. 
Figure  12  shows  the  effect  of  a  thermally  imperfect  interface  on  the  local  work 
lost  for  Qgic  =  1®-  As  interface  'degrades',  i.e.  decreasing  the  heat 
exchange  between  the  slabs  decreases,  and  the  work  lost  decreases  as 
expected. 

Using  Eq.  (108)  the  local  entropy  produced  per  cycle  in  each  layer 
ASjix,Gi),j  =  1,2,  is  given  by 


A  /  \  it  +y„y»,)  p  p  /•lo'Tv 


n=l  ro=l  ' 


As  a  check  on  the  calculations,  the  results  of  Kinra  and  Milligan  (1994)  for 


two  semi-infinite  rods  in  perfect  thermal  contact  were  reproduced  by  letting 
Qgic  Eq.  (127)  become  sufficiently  large  so  that  the  heat  conduction  was 
confined  near  the  interface.  A  normalized  local  entropy  produced  %  is 
introduced  as 
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Normalized  Position  ^ 


Figure  12.  Normalized  work  lost  \)/  as  a  fimction  of  the  normalized  position  ^ 
for  problem  (1)  with  an  imperfect  thermal  interface,  for  the  composite  SiC-Al, 
for  four  normalized  thermal  contact  conductances  for  a  normalized 
frequency  Qgic  for  a  SiC  volume  fraction  of  )^  =  0.5. 
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X=;,  %,  '  (128) 

((acTa)  /c)^ 

A  physical  feel  for  the  denominator  was  given  in  the  paragraph  containing 
Eq.  (126). 

For  the  case  of  a  thermally  perfect  interface,  %  is  plotted  in  Figure  13 
versus  the  normalized  position  ^  with  the  normalized  frequency  Qjic  ^  ^ 
parameter.  Note  that  the  normalization  of  Aw  and  was  done  in  such  a 
manner  that  setting  vj/  =  x  corresponds  to  Aw=7^As,  i.e.,  one  unit  of 
normalized  entropy  produced  corresponds  to  one  unit  of  normalized  work  lost. 
Accordingly,  i|/  (Figures  11  and  12)  and  x  (Figure  13)  are  plotted  on  an 
identical  scale.  As  in  Section 3.1.2,  pointwise  or  Aw^tT^As.  As 

expected  the  entropy  produced  is  maximum  at  the  interface  since  the 
temperature  gradients  are  greatest  there  (see  Figure  8  and  Eq.  (11)),  and  is 
zero  at  the  adiabatic  planes,  ^  =  0  and  ^  =  1.  Figure  14  shows  the  effect  of  a 
thermally  imperfect  interface  on  the  local  entropy  produced  for  Qsic  ~  ^ 

the  interface  'degrades',  i.e.  ^  decreases,  the  heat  exchange  between  the  slabs 
decreases,  and  the  entropy  production  decreases. 

The  maximum  elastic  stored  energy  W  per  unit  volume  during  a  cycle  is 
given  by 


(129) 
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adiaiatic  plane  adiabadc  plane 


Normalized  Position,  ^ 


Figure  13.  Normalized  entropy  produced  x  as  a  function  of  the  normalized 
position  ^  for  problem  (1)  with  a  perfect  thermal  interface,  for  the  composite 
SiC-Al,  for  four  values  of  the  normalized  frequency  Qgicj  for  a  SiC 
volume  firaction  oiVj  =  0.5. 
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adiabatic  plane 


Interface 


adiabatic  plane 


0.00  0.25  0.50  0.75  1.00 


Normalized  Position,  ^ 


Figure  14.  Normalized  entropy  produced  xas  a  function  of  the  normalized 
position  4  for  problem  (1)  with  an  imperfect  thermal  interface,  for  the 
composite  SiC-Al,  for  four  normalized  thermal  contact  conductances  for  a 
normalized  frequency  Qgic  =  for  a  SiC  volume  fraction  of  =  0.5. 
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where  v  is  Poisson's  ratio.  Substituting  Eqs.  (103)  and  (129)  into  Eq.  (4)  and 
using  the  results  of  this  section,  the  specific  damping  capacity  'P  averaged 
over  the  volume  of  the  unit  cell  for  both  problems  (1)  and  (2)  is 


'P  =  'P„ 


n=l  +7/1 


a]  E, 


(130) 


where  Q,  = /^(/la/i  +  ^Faj/a,).  A  modulus  of  elastothermodynamic 
damping  has  been  defined  as  =  2nEa^T„IC.  This  parameter  first  appeared 
in  Zener's  work  (1937)  when  he  analyzed  the  thermoelastic  damping  of  a 
homogeneous  beam  imdergoing  flexural  vibrations.  Zener  called  it  the 
relaxation  strength.  The  damping  results  of  Milligan  and  Kinra  (1993)  for  a 
linear  inclusion  in  an  infinite  matrix  were  reproduced  using  the  above 
equation  by  letting  the  parameter  ajb  become  sufficiently  small  for  a  fixed 
value  of  Q,. 

The  damping  of  the  SiC-Al  composite  with  a  thermally  perfect  interface 
is  plotted  in  Figure  15  as  a  fimction  of  Qgic  for  various  SiC  volume  Auctions, 
Vj .  The  equilibrium  temperature  7^  has  been  taken  to  be  300  K.  As  Qgic  ® 
essentially  isothermal  conditions  are  obtained  and  4^-»0.  At  the  other 
extreme,  as  Qgic  -» <»  essentially  adiabatic  conditions  are  obtained  and,  once 
again,  T  ->  0.  In  between  these  two  extremes,  the  peak  damping  occurs  at  a 
frequency  which  depends  on  the  volume  fraction.  The  maximum  damping  is 
slightly  more  than  0.5%  at  « 10  and  «  0.5.  Now,  for  a  homogeneous 
medium  subjected  to  a  homogeneous  stress  field,  the  damping  is  zero. 
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logio  QsiC 

Figure  15.  Specific  damping  capacity  'P  as  a  fimction  of  the  normalized 
frequency  Q^c  for  various  SiC  volume  fractions  Vj,  for  the  composite  SiC- 
Al,  for  problem  (1)  with  a  thermally  perfect  interface. 
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Accordingly,  in  Figure  15  the  damping  goes  to  zero  for  ->  0  and  ^  1  at 
any  frequency.  The  effect  of  a  thermally  imperfect  interface  on  the  specific 
damping  capacity  is  shown  in  Figure  16  for  Qsic=10  ^nd  Vy=0.5.  As 
expected  from  the  previous  discussion,  as  ^  ^  0  (thermally  isolated  slabs), 
'F  ^  0;  'F  increases  monotonically  with  as  ^  oo,  SF  increases  to  its  value 
for  a  thermally  perfect  interface. 

Finally,  using  the  Simplex  method  of  multivariable  optimization  (Press,  et 
al,  1989,  pp.  289-293),  the  optimum  damping  values  for  various  inclusion 
materials  in  an  Al,  Mg,  and  Ti  matrix,  varying  Vj-  and  Q,,  were  calculated  for 
bofii  problems  (1)  and  (2).  These  results  are  given  in  Table  II.  The 
thermomechanical  properties  of  the  materials  are  given  in  Table  I. 


Table  2  Optimum  specific  damping  capacity  for  problems  (1)  and  (2)  with  a 
thermally  perfect  interface  for  various  inclusion  and  matrix  materials,  varying 
inclusion  volume  fraction  Vj  and  the  normalized  frequency  fij. 


Problem  (1) 

Problem  (2) 

Vf 

Qj 

max 

Vf 

Qj 

'F 

max 

SiC-Al 

0.51 

5.1 

0.0052 

0.26 

1.9 

0.0002 

AI2O3-AI 

0.54 

5.3 

0.0027 

021 

2.8 

0.0030 

SiC-Mg 

0.50 

5.5 

0.0074 

0.22 

1.4 

0.0001 

AljOj-Mg 

0.51 

5.7 

0.0005 

0.22 

2.3 

0.0003 

SiC-Ti 

0.57 

3.2 

0.0008 

0.34 

0.9 

0.0000 

Al,03-Ti 

0.54 

3.8 

0.0001 

031 

2.0 

0.0008 
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Figure  16.  Normalized  specific  damping  capacity  as  a  function  of 

the  normalized  film  coefiScient  for  problem  (1)  with  a  thermally  imperfect 
interface,  for  the  composite  SiC-Al,  for  a  normalized  frequency  Qgjc  =  10,  and 
for  a  SiC  volume  fraction  of  1^  =  0.5. 


3.4.2  Elastothermodynamic  Damping  in  a  Symmetric  Three-Layer  Plate 
Consider  a  symmetric  three-layer  plate  with  thermally  perfect  interfaces 
as  shown  in  Figure  17.  In  the  following,  all  externally  applied  loads  are 
assumed  to  be  time-harmonic  and  in  phase  with  each  other,  and  the  factor  e'“' 
is  implied,  as  is  customary.  It  is  assumed  that  the  deformations  in  the  plate  are 
adequately  described  by  the  classical  (thin)  laminate  theory  based  on  the  well- 
known  Kirchoff  hypothesis  (CLT).  Since  the  plate  is  symmetric  in  both  the 
material  properties  and  the  geometry,  there  is  no  coupling  between  bending 
and  extension.  Herein,  only  bending  is  considered.  Then,  given  the 
bending/twisting  moments  M^{x,y),  My{x,y),  and  M^{x,y),  the  curvatures 
K,(x,y),  Kj,(x,y),  and  K^{x,y),  may  be  obtained  by  using  the  CLT  (see 
Jones,  1975,  for  example).  Our  sign  convention  is  that  positive  moments 
produce  positive  curvatures  (see  Figure  17);  it  is  different  from  that  of  Jones 
(1975).  With  CTjj  =  =  CT^  =  0,  the  nonzero  stresses  are  given  by 

cy„=-j^(K,  +  VKjz,  a^=-j^(K^+vK,)z, 

(131) 

so  that 

=-Y^(k,+kJz 


(132) 
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Note  that  the  hydrostatic  stress  field,  <5^,  is  independent  of  the  twist 
curvature,  ,  and,  therefore,  of  the  twisting  moment,  .  In  keeping  with 
the  assumptions  in  the  CLT,  the  gradients  of  temperature  in  the  x  and  y 
directions  are  very  small  compared  to  the  gradients  in  the  z-direction.  From 
Eq.  (11),  the  production  of  entropy  depends  upon  the  square  of  these 
temperature  gradients.  Therefore,  for  the  purpose  of  calculating  the  total 
damping  of  the  laminate,  heat  conduction  in  the  x  and  y  directions  is 
neglected. 

Note  that  in  Eq.  (132)  is  an  odd  fimction  of  z.  Since  the  operator  in 
Eq.  (15)  is  symmetric  in  z,  it  follows  that  F"  is  an  odd  fimction  of  z,  and  V  =  0 
at  z  =  0.  Thus,  it  is  sufficient  to  consider  only  the  upper  half  of  the  plate, 
z>0.  In  the  following  the  subscript  (  ),  will  refer  to  the  inner  layer  of 
thickness  2a,  and  the  subscript  (  )2  will  refer  to  the  two  outer  layers  of 
thickness  {b-a);  the  total  thickness  of  the  plate  is  2b  (see  Figure  17).  The 
volume  firaction  of  the  outer  layers  is  Vj-  =  l-alb.  Using  the  method  of 
Ozisik  (1980,  pp.  294-334),  the  eigenfimctions  for  each  of  the  laminae  are 

K(z)  =  sm{y„zla)  (133) 

<t>2»  (^)  =  A,  <^os{r\„zlb)  +  B„  sin(ii„z/b)  (134) 

where  A„  and  are 


A=- 


iCcosY„  cosr|„ 
sin(Ti„(f-l))  ’ 


i:cosY„sinr| 

sin(Ti„(f-l)) 


(135) 


and  y„,  t|„,  and  K  were  defined  in  Eq.  (1 12).  The  transcendental  equation  for 
the  eigenvalues  is 
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-siny„ 

cos(Ti„a[/6) 

sin(n„a/^)) 

/rcosy„ 

sin(ii„a/Z>) 

-cos(ti„a/6] 

1=0 

(136) 

0 

-sinTi„ 

cos-n„ 

Using  exactly  the  same  procedure  as  in  Section  3.4.1,  the  specific 
damping  capacity  for  the  three-layer  composite  plate  was  foimd  to  be 
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— (cos(ti«  )  -  cos(Ti„a/6))  -I-  sin(ri„ )  -  ^  sin(ri„a/6) 

+B„ 

— (sin(Ti„)  -  sin(Ti„a/6))  -  cos(ti„  ) +^cos(r\„a/b) 

r  1  1  ^ 

A  =r-^sin(2y„) 

2  4y„ 


-  [cos(2ii„ )  -  cos(2Ti„a/6)] 


(140) 


The  damping  of  the  zinc-steel-zinc  composite  plate  is  plotted  in  Figure  1 8 
as  a  function  of  for  various  values  of  Vj  (the  volume  fraction  of  zinc), 
with  the  cinwatures  =  0,  k,  =  =  constant,  and  7^  =  300iC.  For  a  given 

value  ofVy,as  —> 0  essentially  isothermal  conditions  are  obtained,  and 
'F->0.  At  the  other  extreme,  as  away  from  the  interfaces, 

essentially  adiabatic  conditions  are  obtained,  and,  once  again  'F  ^  0.  For 
comparison,  T  for  an  uncoated  steel  plate  [v^  =0.o)  is  also  shown;  the 
maximum  damping  is  =0.025.  At  a  "very  small  cost"  of  adding  20% 
zinc  (by  volume),  'F^a,^  increases  to  0.055,  i.e.  by  about  100%.  When 
Vf  =0.4,  'F^ax  =  0-l»  i-C-j  an  increase  of  about  400%. 
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Figure  18.  Specific  damping  capacity  'F  versus  the  normalized  frequency 
for  the  symmetric  three-layer  plate  zinc-steel-zinc  with  thermally  perfect 
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The  special  case  of  a  homogeneous  plate  may  be  obtained  from  Eq.  (137) 
by  letting  Vj  ->  0.  After  some  careful  analysis,  the  following  result  was 
obtained 


'F  = 


1+v/3(k) 


Q. 


(141) 


1  -  V  /,  (K)  L  ■ "  (2n  + 1)^  +  (In  + 1)^ 

where  an  unsubscripted  normalized  frequency  has  been  defined  as 
Q.  =  (i)C{2a)^  jn^k,  so  that  the  result  will  be  directly  comparable  to  a  classic 
result  by  Zener  (1937).  Note  that  2a  is  the  thickness  of  the  homogeneous 
plate.  The  expression  in  the  square  brackets  in  Eq.  (141)  is  precisely  the 
damping  of  a  homogeneous  Euler-Bemoulli  beam  subjected  to  in-plane 
bending  (Zener,  1937).  Thus,  the  following  elegantly  simple  result  follows 
immediately.  For  a  homogeneous  plate. 


(142) 

where  X,  =  (1+v)/3(k)/(1-v)/,(k).  Consider  the  special  case,  k^=0, 
K,  =  =  constant,  and  v  =  1/3;  then  X  =  3.  Therefore,  'Fp,j,e  =  i.e., 

the  damping  in  a  plate  subjected  to  equal  biaxial  curvatures  is  three  times  the 
damping  in  a  beam  subjected  to  the  same  curvature  in  the  plane  of  the  beam. 
Furthermore,  if  Kj^=k^  =  0  then  A,  =  2,  i.e.,  when  subjected  to  identical 
ctirvatures,  the  plate  damping  is  twice  the  beam  damping.  Zener  (1937)  made 
a  clever  observation  that  the  series  in  Eq.  (141)  is  accurately  represented  by 
merely  its  first  term  (n  =  0),  and  that  96/71^  =0.9855  =  1.  Then,  Eq.  (141) 
reduces  to  the  well-known  form  for  a  standard  anelastic  solid 


Q 

1+Q^ 


(143) 


3.43  Elastothermodynamic  Damping  in  a  Hollow  Two-Layer  Concentric 
Sphere 

Consider  a  random  particulate  composite  reinforced  by  hollow  spheres  of 
outer  radius  a  with  the  volume  fraction  of  inclusions  given  by  1^ .  As  a  first 
approximation  to  the  elastothermodynamic  damping  of  such  a  composite,  we 
calculate  the  damping  in  a  hollow  two-layer  concentric  sphere  with  the 
inclusion  occupying  the  space  h<r<a  and  the  matrix  occupying  the  space 
a<r<b  such  that  {ajb^  =Vj  as  shown  in  Figure  19.  Thus,  the  volume 
fraction  of  the  inclusion  in  the  two-layer  composite  sphere  is  equal  to  the 
volume  fraction  of  inclusions  in  the  random  particulate  composite.  The 
surfaces  r  =  h  and  r  =  6  are  assumed  to  be  adiabatic  and  the  interface  r  =  a 
may  be  thermally  imperfect.  The  outer  surface  r  =  b  is  subjected  to  a  time- 
harmonic  radial  stress  =  ae'“' .  Unless  otherwise  stated,  the  subscripts  (  ), 
and  (  )2  will  refer  to  the  inclusion  and  matrix,  respectively. 

It  is  assumed  that  the  frequency  of  vibration  is  sufficiently  small  so  that 
the  inertia  effects  may  be  neglected  and  the  dynamic  stress  field  may  be 
approximated  by  the  corresponding  static  stress  field  times  The  static 
elasticity  problem  may  be  solved  using  standard  methods  (Little,  1973). 
Interestingly,  the  hydrostatic  stress  is  constant  in  the  inclusion  and  the  matrix. 
The  pertinent  results  are 


•IIU 


with  dimensions  shown. 


— jw  sphere 
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+  6^,D^  (Vo)’)(l  -  (hla)%lbf  + 

+(f  +6lx,D^{a/bf\l-{‘‘/tf) 

where  K  is  the  bulk  modulus,  is  the  shear  modulus,  and  the  constants 
D, ,  ,  Dj ,  D4  are  given  by 


°  I  uia -i!iL((„/*)> -1)],  D, 

3K,  '  ’  ’  ’\  ^  detl^  3^J 


A  =  A  +0).  A  =  A  (i-tS) 


where 


det  =  -  \+^{alhy  A\xJ^{alhf l  +  ^(a/b)^  4^,((a//i)^-l)  (147) 

The  eigenfunctions  for  each  layer  may  be  determined  using  the  method  of 
Ozisik  (1980,  pp.  294-334).  The  eigenfunction  for  each  region  is 

‘l>i„  -  (lA)sin(y„x/fl)  +  5,„(l/x)cos(y„x/a)  (148) 

<t>2«  =  4„(lA)sin(Tl„x/6)+52„(l/x)cos(Ti„x/6)  (149) 

where  the  constants  ,  A2„ ,  B2„  are  given  by 

^l„=  — .  4»=V"44(«12«21-«22«1i)»  ■B2«  =  ^ («12%  " «22«1  1 )  O^O) 

a,2  A  A 

where  A  =  ai2(a23a44 -a43n24),  and  y„,  ri„  are  given  by  Eq.  (112).  The 


transcendental  equation  for  the  eigenvalues  p„  is 
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«i2  0  0 


^21  ^22  ^23  ^24 

aji  ^33  ^34 

0  0  a43  a44 


(151) 


where  the  constants  ,  i,j  e  (1,2, 3, 4}  are  given  by 
an={yr,h/a)cos{y„hla)-sm{y„h/a),  0,3  = -(y„/z/a)sin(y„V«)-cos(Y„V«) 
«2i  =feA2)(Y«cosY„-sinYj,  a22  =-(^iA2)(Y»sinY„+cosYj 

«23  =  -(T1„«/^)cos(Ti„a/Z>) + sin(Ti„a/Z>) 

«24  =  {‘nn^/b)  sm{T\„a/b) + cos{T]„a/b) 

«3i  =(^iA2)(-Y„cosy„  +sinY„)-(/a/^2)sinY„ 

«32  =(^iA2)(Y„sinY„  +cosYj-(/aA2)cosY„ 

^33  =  {fa/k2)sm{i]„a/b),  =  {falk^)c,oir\„alb) 

«43  =  'Hn  cos(TiJ-sin(TiJ,  a44  =  -n„  sin(Ti„)-cos(TiJ  (152) 
Using  exactly  the  same  procedure  as  in  Section  3.4.1  the  specific  damping 
capacity  for  the  entire  spherical  composite  was  found  to  be 


«  r\  A, 2 


_ ^ 

]:{alb)\\-{hla)^)  +  ^^{DjD^f{alb)\hlaf[\-{hla)^) 

1  Z  J  Aj 


(153) 


76 


where  Q,  =(oC\a^  fk^,  a  modulus  of  damping  for  this  problem  is  defined  as 
=  TzTgKa^  fC,  and  Q„  is  defined  as 


with 


0,= 


(154) 
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h  “  4"J[sin(2Ti„)-sin(2Ti„a/6)] 

-  ^  AnK  [cos(2ti„  )  -  cos{2r\„a/b)] 

^  *1/1 

(155) 

The  damping  is  plotted  in  Figure  20  for  the  composite  AljOj-Al  with  a 
thermally  perfect  interface  for  various  inclusion  volume  fractions  with  hja  =  0 
(solid  sphere).  Even  though  the  theoretical  maximum  voliime  fraction  of  an 
actual  particulate  composite  containing  spherical  inclusions  of  uniform  size  is 
0.74,  the  damping  is  plotted  for  >0.7  for  completeness.  The  maximum 
damping  is  approximately  0.01  for  «0.5  and  Qj  »30  (cf.  Section  3.4.1, 
Table  2). 

Figure  21  shows  the  effect  of  the  cavity  ratio  hja  on  the  damping  of  the 
AljOj-Al  composite  with  a  thermally  perfect  interface  for  Fy=  0.5  and 
several  normalized  frequencies.  The  damping  of  the  composite  is  relatively 
unaffected  by  the  hollow  interior  within  the  range  0  <  hja  <  0.5.  The  limiting 
case  hla  =  l  corresponds  to  a  spherical  shell.  For  this  case,  since  the 
hydrostatic  stress  is  constant,  tiie  thermoelastic  effect  is  homogeneous,  and 
hence  there  is  no  conduction  of  heat.  Therefore,  the  damping  is  zero. 
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Figure  20.  Specific  damping  capacity  'F  as  a  function  of  the  normalized 
frequency  Q,  for  the  spherical  composite  AljOj-Al  with  a  thermally  perfect 
interface,  for  various  AI2O3  volume  fractions  Vj ,  with  hfa  =  0. 


Specific  Damping  Capacity, 


79 


cavity  ratio,  h/a 


Figure  21.  Specific  damping  capacity  T  as  a  function  of  the  cavity  ratio  hja 
for  the  spherical  composite  AljOj-Al  with  a  thermally  perfect  interface  for 
various  normalized  frequencies  Q,  and  =0.5. 
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Finally,  using  the  Simplex  method  of  multivariable  optimization  (Press,  et 
al,  1989,  pp.  289-293),  the  optimum  damping  values  for  SiC  and  AljOj 
inclusions  in  an  Al,  Mg,  and  Ti  matrix,  varying  Vj  and  Qj  with  hja  =  0,  were 
calculated.  These  results  are  given  in  Table  3.  The  thermomechanical 
properties  of  the  materials  are  given  in  Table  1. 

Table  3.  Optimum  specific  damping  capacity  for  the  composite  sphere  with  a 
thermally  perfect  interface  and  hja^Q,  for  various  inclusion  and  matrix 
materials,  varying  volume  fraction  Vj  and  the  normalized  frequency  Q,. 


Vf 

Q, 

'P 

^  max 

SiC-Al 

0.44 

36 

0.028 

AI2O3-AI 

0.47 

32 

0.011 

SiC-Mg 

0.46 

46 

0.042 

AljOj-Mg 

0.47 

46 

0.022 

SiC-Ti 

0.44 

15 

0.004 

AljOj-Ti 

0.44 

21 

0.000 
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4.  ELASTOTHERMODYNAMIC  DAMPING  IN 
CRACKED  MEDIA 

In  this  section,  a  Griffith  crack  subjected  to  a  time-harmonic  loading  in 
Modes  I,  II,  and  III  is  considered.  An  approximate  analysis  is  given  for  the 
temperature  field,  the  work  lost,  and  the  entropy  produced  in  the  vicinity  of 
the  Griffith  crack. 

4.1  STRESS  FIELD 

Consider  an  infinite  plate  of  an  isotropic  homogeneous  material  (in  either 
plane  stress  or  plane  strain)  containing  a  Griffith  crack  occupying  the  line 
-a<x<a,  y  =  Q.  It  is  assumed  that  the  fi-equency  of  vibration  is  sufficiently 
small  (or,  alternatively,  the  wavelength  is  very  large  compared  to  the  crack 
length)  so  that  the  inertia  effects  may  be  neglected  and  the  dynamic  stress 
field  may  be  approximated  by  the  corresponding  static  stress  field  times 
The  hydrostatic  component  of  the  static  stress  field  is  given  by  (Westergaard, 
1939;  Paris  and  Sih,  1965) 


Mode  I: 

z 

<3 u  =  2K(y  Re  -  = 

(156) 

Mode  II: 

a  It  =2KTlm  ■ , . 

**  ^2 

Vz  —a 

(157) 

Mode  III: 

q 

II 

O 

(158) 
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where  z  =  x+iy,  k  =  1  for  plane  stress  and  1  +  v  for  plane  strain,  v  is 
Poisson's  ratio,  and  cj  and  t  characterize  the  stress  field  at  infinity: 
Cii  =  <J22  =  C7  and  a, 2  =  0  for  Mode  I;  c„  =  022  =  0  and  <3^^  —  for  Mode  II. 

4.2  MODE  III 

First,  the  trivial  case  of  Mode  III,  where  is  zero  everywhere,  is 
dismissed.  Since  is  zero  everywhere,  the  thermoelastic  effect,  Eq.  (1), 
disappears.  There  are  no  changes  in  temperature,  no  entropy  produced  and  no 
work  converted  into  heat.  Attention  is  now  turned  to  Modes  I  and  II.  Since 
closed-form  solutions  were  obtained  for  Mode  II  but  not  for  Mode  I,  Mode  II 
results  will  be  given  first. 

4.3  MODE  II 

4.3.1  Temperature  Field 

The  governing  equation  for  the  fluctuating  temperature  field  is  given  by 
Eq.  (15), 

V'r-^F  +  ig(x)  =  0  (159) 

where  the  inhomogeneous  stress  term  has  been  written  as  a  heat  generation 
term  g(x)  given  by 

g{x)  =  -(/© )  aT„  cTtt  (x)  (160) 

where  <3 (x)  is  given  by  Eq.  (157).  In  Mode  II  (x)  is  an  odd  function  of 
X  and  y.  Therefore,  with  the  exception  of  the  crack  faces,  -a  <  x  <  a,  >»  =  0*, 
Cjyt  =0  on  the  planes  x  =  0  and  3^  =  0.  Since  the  operator  in  Eq.  (159)  is 
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symmetric  in  x  and  y,  F  is  an  odd  function  of  x  and  y  and,  with  the  exception 
of  the  crack  faces,  -a<x<a,y  =  0*,V  =  Qonx  =  0  and  y  =  0.  The  crack 
faces  are  taken  to  be  adiabatic,  i.e.  dVfdy  =  0.  Therefore,  instead  of  the  full 
plane  shown  in  Figure  22(a)  it  is  sufficient  to  consider  the  quarter  plane  shown 
in  Figure  22(b)  with  the  following  thermal  boundary  conditions 

(1)  F  =  0  x  =  0,0<y<oo 

(2)  F  =  0  y  =  0,a<x<«>  (161) 

(3)  —  =  0  y  =  0^0<x<a 
dy 

Attempts  to  obtain  an  exact  solution  to  this  problem  led  to  intractable 
mathematics.  An  approximate  solution  is  now  developed.  Accordingly,  a 
number  of  assumptions  intended  to  reduce  the  algebraic  effort  will  be 
introduced.  Once  the  approximate  solution  is  obtained,  the  reasonableness  of 
the  assumptions  introduced  en  route  will  be  demonstrated. 

The  temperature  field,  Eq.  (159),  is  driven  by  via  the  heat  generation 
term,  Eq.  (160).  The  rate  of  entropy  production  is  driven  by  the  square  of  the 
spatial  gradients  of  temperature,  Eq.  (11).  Sufficiently  far  away  from  the 
crack  0,  VT  0,  and  the  elastothermodynamic  effects  (work  lost  and 
entropy  produced)  disappear.  Accordingly,  attention  is  focused  on  the  region 
surrounding  the  crack.  A  polar  coordinate  system  (r,0)  is  introduced  at  the 
right  crack  tip  z  =  a  as  shown  in  Figure  23.  In  Figure  24  the  exact 
contours  are  plotted  for  the  region  0  <  x  <  2a,  0  <  y  <  a,  and  in  Figure  25  the 
exact  is  plotted  as  solid  lines  for  0<r<a  and  0  =  0,7t/4,7i/2,37i/4,7t. 


Figure  22.  (a)  A  Griffith  crack  of  length  2a  in  Mode  II  loading  with  a 
rectangular  coordinate  system  (x,y).  The  crack  faces  are  assumed  to  be 
adiabatic.  Because  Cy,  is  an  odd  function  of  x  and  y,  V  =  0  on  x  =  0,  and  on 
^  =  0  except  on  the  crack  faces,  (b)  The  reduced  quarter-plane  heat- 
conduction  problem. 
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For  a  positive  applied  shear  stress,  the  hydrostatic  stress  field  in  the  first 
quadrant  (x  >  0,y  >  0)  is  compressive;  accordingly,  in  Figure  25  negative 
has  been  plotted.  Moreover,  the  Griffith  crack  causes  a  mode  conversion  of 
the  remotely  applied  pure  shear  stress  =  0)  to  a  hydrostatic  stress  field 
near  the  crack.  As  expected  on  the  basis  of  symmetry,  on  the  ray 

0  =  0,  and  =  0  at  the  point  r  =  a,B  =  Tt  (x  =  y  =  0).  By  virtue  of  the  fact 
that  is  an  odd  function  of  x,  for  0  =  7i,  as  r  increases  from  a  to  2a 
(-a  <x<0)c^  increases  from  zero  to  plus-infinity.  (It  is  for  this  reason  that, 
in  the  next  paragraph,  our  Laurent  series  expansion  of  ajj*  will  be  limited  to 
r  < a.)  Thus,  in  Mode  II  the  Griffith  crack  localizes  the  hydrostatic  stress 
field  to  the  region  between  the  crack  tips.  Later  it  "will  be  shown  that  the  exact 
opposite  is  true  for  Mode  I:  the  Griffith  crack  localizes  the  hydrostatic  stress 
field  to  the  region  exterior  to  the  crack  tips. 

Next,  an  approximation  to  is  constructed.  Eq.  (157)  was  expanded 
into  a  Laurent  series  in  z/2a  about  the  tip  z  =  a.  In  linear  elastic  fracture 
mechanics  (LEFM),  a  characterization  of  the  crack-tip  stress  field  is  generally 
considered  adequate.  This  is  not  the  case  here.  Here,  interest  is  in  the 
conversion  of  mechanical  energy  into  heat  throughout  the  plate  due  to  the 
presence  of  a  crack.  A  solution  is  sought  that  is  valid  in  as  large  a  region 
surrounding  the  crack  as  possible.  Accordingly,  one-,  two-,  three-,  and  four- 
term  expansions  were  compared  with  the  exact  Gu^.  The  expansion  was 
limited  to  the  disk  r<a,  -K<Q<n.  It  was  found  that  the  four-term 
expansion  is  not  a  significant  improvement  over  the  three-term  expansion 
(EssentMly,  the  error  changes  sign).  Furthermore,  retaining  each  additional 
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term  entails  a  significant  amount  of  additional  algebraic  effort.  Therefore,  the 
series  was  truncated  at  three  terms.  In  the  disk  0 <r J a  <\,  -n<0  <7:,  the 
exact  given  by  Eq.  (157)  will  be  approximated  by 


=2kt| 


sin(e/2)  + j  sin(e/2)  ^  sin(30/2)f 

2  V  r  4  V2a  16  {2a J 


0<r  <a 


(162) 


-71  <  0  <  71 

This  approximate  is  compared  with  the  exact  in  Figure  25.  The 
comparison  is  considered  a  good  engineering  approximation.  (At  the  risk  of 
stating  the  obvious,  outside  the  disk  r  =  a,  Eq.  (162)  is  a  poor  approximation 
to  Eq.  (157).  For  example,  as  noted  earlier,  on  0  =  71,  a  <r  <2a,  the  exact 
increases  from  zero  to  plus  infinity  whereas  the  approximate 

Ott^(9/8)KT.) 

Next,  approximations  to  the  heat-conduction  problem  shown  in 
Figure  22(b)  are  considered.  A  normalized  frequency  Q  and  a  diffusion  length 
/  are  introduced  as 


Q  =  (i>Cayk,  l  =  ^2kl(oC  (163) 

so  that  Q  =  2{afiy  or  /  =  a(2/Q)’^^ .  A  "physical  feel"  for  the  diffusion  length 
may  be  obtained  as  follows.  Consider  a  semi-infinite  rod  occupying  0<x<<x> 
at  an  equilibrium  temperature  2^,  subjected  to  a  time-harmonic  temperature 
V  =  at  the  end  x  =  0.  In  the  interior  the  temperature  variation  is  given 
by  F(x)  =  F^exp(-(l  +  i)x//).  Thus,  at  x  =  l,  |r/]^|  =  l/e,  and  so  on. 
Returning  to  the  problem  at  hand,  observe  that  the  source  term  g(x)  in 
Eq.  (159)  is  square-root  singular  at  the  crack  tip.  When  I  is  comparable  to  the 
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half-crack  length  a,  Q  =  0{1),  a  strong  thermal  interaction  between  the  crack 
tip  and  the  isothermal  wall  at  ;c  =  0  is  anticipated.  In  order  to  capture  this 
physics  one  must  obtain  the  exact  solution;  the  approximate  solution  will  not 
do  that.  Suppose  Q  is  large  compared  to  one,  say  Q  >  10,  then  (a/l)  >  2.24. 
Therefore,  a  rather  weak  interaction  between  the  crack  tip  and  the  isothermal 
wall  at  X  =  0  is  anticipated.  In  the  following  it  is  assumed  that  there  is  no 
interaction  between  the  crack  tip  and  the  wall  at  x  =  0.  This  allows  for  the 
construction  of  a  much  simpler  heat-conduction  problem  associated  with  the 
actual  heat-conduction  problem  shown  in  Figure  22(b).  Appropriately,  it  will 
be  called  the  "associated  heat-conduction  problem";  this  is  shown  in 
Figure  26.  Since  the  wall  is  "thermally  unaware"  of  the  crack  tip  at  x  =  a,  it 
can  be  moved  from  x  =  0  to  x  =  -co  with  impunity;  concurrently,  the 
adiabatic  line  0  <  x  < a  stretches  out  to  the  adiabatic  line  -oo  <x<a. 

Attention  is  now  turned  to  the  heat  generation  term  g(r,9).  In  the 
interior  of  the  half-disk,  r<a,  O<0<7t,  g(r,9)  for  the  associated  problem  is 
taken  to  be  that  given  by  the  three-term  expansion,  Eqs.  (160)  and  (162). 
Returning  to  Figure  25,  at  r  =  a,  is  rather  small;  more  significantly,  the 
gradients  of  0*^  are  small.  Therefore,  in  the  associated  problem  we  set 
g(r,9)  =  0  in  a<r<oo,  0^9 <n.  Once  the  solution  to  the  associated 
problem  is  obtained,  the  solution  in  the  left  quarter-plane,  x  <  0,  y  >  0  will  be 
immediately  discarded.  The  solution  in  the  right  quarter  plane,  x  >  0,  y  >  0, 
will  be  taken  to  be  the  solution  for  the  actual  problem  shown  in  Figure  22(b). 
As  noted  earlier,  the  solution  in  the  ftill  plane.  Figure  22(a),  can  be  readily 
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Figure  26.  The  "associated  heat  conduction  problem"  with  zero  heat 
generation  outside  the  disk  of  radius  a.  Mode  11. 
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constructed  using  the  symmetries  in  <y^  and  V.  Finally,  it  is  intuitively 
obvious  that  as  Q  becomes  very  large,  i.e.  as  the  diffusion  length  becomes 
very  small  compared  to  the  crack  length,  the  approximate  solution  will 
become  increasingly  more  accurate. 

The  associated  heat-conduction  problem  is  given  by 

^+l^+l^-£mV  +  jg(r,e)  =  0,  0<e<7i,0<r«»  (164) 

or  r  or  r  59  k  k 


V  =  0  at  0  =  0,  0<r<oo 
dV 


56 


=  0  at  0  =  11,  0<r<oo 


(165) 


where 


^(^e)= 


r-(/co)a7^ a*t(r,6)  O<0<7i,  0<r<a 


(166) 


0  0<6<7t,  a<r<co 

Next,  an  exact  solution  to  the  associated  heat-conduction  problem  is  obtained. 
An  integral  transform-inversion  formula  is  first  developed  to  remove  the  6- 
variable  in  Eq.  (164).  It  is  assumed  that  the  fluctuating  temperature  field 
V(r,d)  may  be  represented  in  the  form 


rfr.e)  =  2;c.«sm((n+l/2)6)  (167) 

where  the  sin((n -i- 1/2)0)  are  the  eigenfunctions  arising  from  solving  the 
corresponding  homogeneous  heat  conduction  problem.  Multiplying  Eq.  (167) 
by  sin((/w+ 1/2)0)  and  integrating  from  0  to  it  with  respect  to  6  yields 

C„(i-)=-F(r) 

It 


(168) 
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where  the  integral  transform  of  V  is  defined  as 

V(r)  =  f"  V(r,e)  sin((«  + 1/2)0)  dQ  (169) 

Jo 

The  inversion  formula  is  obtained  by  substituting  Eq.  (168)  into  Eq.  (167) 
yielding 

r(r,e)  =  -|;F(r)sm((»+ 1/2)0)  (170) 

^  n=0 

Next,  taking  the  integral  transform  of  Eq.  (164),  integrating  the  0- 
derivative  term  by  parts,  and  using  the  boundary  conditions  Eq.  (165)  results 
in 


d^V  ^  1  dV 
dr^  r  dr 


r.  r 


C.  («  +  l/2) 


2  A 


J 


V  =  -Tg^^^>  0<r<«>  (171) 

k 


where  the  transform  of  the  heat  generation  term,  g{r),  is  given  by 


-(za))a7^Kx 


7cV2 


-1 


na  4 
0 


'na 


Q<r<a 
a  <r  <co 


(172) 


where  5„o=0  for  n^O,  5„o=l  for  «  =  0,  and  similarly  6„,  =0  for  n^l, 
6„j  =  1  for  M  =  1.  Eq.  (171)  is  an  inhomogeneous  modified  Bessel's  equation 
of  order  («+l/2).  Two  independent  solutions  to  the  homogeneous  equation 
are  I^„_^y2y{^fiQrfa))  and  where  /(„+i/2)  is  the  modified 

Bessel  function  of  the  first  kind  of  order  (w  +  1/2),  and  ^(„+i/2)  is  the  modified 
Bessel  function  of  the  second  kind  of  order  (n  +  1/2)  (McLachlan,  1955); 
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recall  Q  =  ooa^/Z).  Equation  Eq.  (171)  may  be  solved  using  the  method  of 
Green's  function  (Arfken,  1985)  or  equivalently,  the  method  of  variation  of 
parameters  (Boyce  and  DiPrima,  1986)  resulting  in 

«  (173) 

Substituting  Eq.  (172)  into  Eq.  (173)  and  the  result  into  Eq.  (170)  gives 


for  0<r<a,  O<0£7i.  (174) 

where  a  normalized  fluctuating  temperature  ©  for  Mode  II  has  been 
introduced  as  ©  =  F/(xK7^a/C) .  (A  physical  feel  for  the  denominator  may  be 
developed  as  follows.  If  a  rod  is  subjected  to  a  uniaxial  tension  of  magnitude 
T  imder  adiabatic  conditions,  it  will  cool  by  an  amount  xT^a/C.)  The  Bessel 
functions  appearing  in  Eq.  (174)  are  of  half-integer  order  and  may  be 
expressed  in  terms  of  elementary  functions  as  follows  (Abramowitz  and 


Stegvm,  1970) 
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hi  (^)  =  J—  sinh(z) ,  (z) 


I  TIZ 


V2' 


- 

\2z 


exp(-z) 


= 

V  itz 


1 


cosh(z) — sinh(z) 


K,^(z) 


(175) 


The  integrals  appearing  in  Eq.  (174)  were  calculated  explicitly  and  are  given 
by 


jg^''  /i/2  («  Vifi  )(-• \/«  +  j  ]du=^ 


1  +  - 1  j  cosh(  r/a)  -  sinh(  V^r/a) 


-jL=r— cosh(V/Qr/a)+  f— |  +—  sinhf-v/iQr/a) 
VzT2  a  Kay  zQ  I 


(176) 
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Next,  the  crack-tip  temperature  field  is  examined.  For  a  finite  Q  (finite  t), 
expanding  Eq.  (174)  about  r  =  0  and  retaining  only  the  first  term  yields 


0  =  A(n)^sm(e/2),  r«l 

where  the  complex  function  y4(f2)  is  given  by 


(177) 


(178) 


The  two  bracketed  terms  correspond,  respectively,  to  the  first  (l/Vr)and 
second  (Vr)  terms  of  in  Eq.  (162);  the  third  term  does  not 

contribute  to  A(q).  Moreover,  a  careful  examination  of  the  complete  series 
expansion  of  revealed  that  the  contribution  to  A(q)  of  all  higher-order 
terms  and  higher)  is  also  zero.  Now  suppose  for  a  moment  that  the 
thermal  conductivity  is  zero,  then  from  Eq.  (1),  Q  =  (-<y^^/KT),  i.e.,  together 
with  CTjit  the  temperature  is  also  square-root  singular  at  the  crack  tip.  In  this 
context,  it  is  a  remarkable  observation  that  for  any  ^  >  0,  however  small,  the 
temperature  at  the  crack  tip  is  zero  at  all  frequencies,  Eq.(177). 

Next,  the  high-frequency  behavior  of  the  crack-tip  temperature  field  is 
considered.  Intuitively,  as  Q  becomes  very  large  (/  «a),  it  is  expected  that 
the  elastothermodynamic  effects  will  become  increasingly  localized  near  the 
crack  tip.  It  is  expected  that  the  crack-tip  solution  will  become  independent  of 
a.  As  Q ->  oo,  the  second  bracketed  term  in  Eq.  (178),  which  corresponds  to 
the  second  (^/r)  term  of  Eq.  (162),  vanishes,  and  A{Q)->yl2iCl. 


Equation  Eq.  (177)  reduces  to 
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sin(0/2) ,  r  «l« a  (179) 

It  is  important  to  observe  that  ©  in  Eq.  (179)  corresponds  only  to  the  square- 
root  singular  stress  field  in  Eq.  (162).  In  LEFM  it  is  customary  to  work  with 
dimensional  quantities,  e.g.,  K^.  Accordingly,  Eq.  (179)  is  written  in  a 
dimensional  form  as 

V(r,Q)  =  T,e’^>*--^4^K„J^sm{ej2),  r«l«a  (180) 

where  Kjj  =  Note  that  V  ~  V®  .  In  LEFM,  for  a  fixed  Kji,  the  crack- 

tip  stress  field  for  all  Mode  II  problems  is  identical,  and  is  given  by  the  first 
term  in  Eq.  (162):  =(-2K)sin(6/2)/C;;/(27ir)’^^ .  Similarly,  for  a  fixed 

Kj,  and  ®,  the  high-fi’equency  crack-tip  temperature  field  for  all  Mode  II 
problems  is  identical,  and  is  given  by  Eq.  (180)  (Obviously,  for  plates  made 
up  of  different  materials,  the  temperature  field  depends  on  the  thermal 
property,  a/V^.). 

In  polar  form  ©  =  ||©||exp(i  O),  where  ||©||  is  the  magnitude  and  <I)  is  the 
phase  of  ©.  Returning  to  Eq.  (174),  in  Figure  27  the  magnitude  of  ©  is  plotted 
for  0<r  <a,  0  =  O,7t/4,7c/2,37i/4,7t,  and  =  1,10,100,1000.  It  is  recognized 
that  the  approximate  solution  may  begin  to  lose  accuracy  below  Q  =  10. 
Nevertheless,  ©  is  plotted  for  =  1  to  give  a  "feel"  for  it.  The  fluctuating 
temperature  field  obtained  using  the  two-term  expansion  for  is  also 
included  for  comparison.  Physically,  at  low  frequencies  heat  has  sufficient 
time  to  diffuse  during  each  cycle  so  that  essentially  isothermal  conditions  are 


Magnitude  of  the  Normalized  Magnitude  of  the  Normalized 

Fluctuating  Tenperature,  ||@||  Fluctuating  Tenperature,  ||©|| 
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Normalized  Radial  Coordinate,  r/a  Normalized  Radial  Coordinate,  r/a 
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Figure  27.  Magnitude  of  the  normalized  fluctuating  temperature  as  a  function 
of  the  normalized  radial  coordinate  for  9  =  0, 7t/4,7i/2,37c/4,7t .  (a)  Q  =  1;  (b) 
Q  =  10;  (c)  =  100;  (d)  Q  =  1000.  The  fluctuating  temperature  at  the  crack 

tip  is  always  zero.  Mode  II. 
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obtained.  This  case  is  evidenced  by  Q  =  1  in  Figure  27.  Conversely,  at  high 
frequencies,  heat  has  very  little  time  to  diffuse  and,  away  from  the  crack  tip, 
essentially  adiabatic  conditions  areobtained.  Under  adiabatic  conditions, 
from(l),  ©  =  (-ajj/KT).  This  explains  the  quantitative  agreement,  away 
from  the  crack  tip,  between  0  in  Figure  27  and  /kt)  in  Figure  25, 
particularly  at  =  1000.  Close  to  the  crack  tip,  the  temperature  gradients  in 
the  radial  direction  are  square-root  singular  (see  Eq.  (177))  and,  therefore, 
truly  adiabatic  conditions  are  never  attained  even  at  very  high  frequencies. 
This  explains  the  difference  between  0  and  (-a*;i/KT)  close  to  the  crack  tip. 
This  argument  is  now  explored  mathematically.  For  any  r  e(0,fl),  it  is  easy 
to  show  from  Eqs..(174)  and  (176)  that  lim0(r,0)  =  (-CT*ji/KT).  However, 
as  r->0,  becomes  unbounded,  whereas  0  is  a  bounded  continuous 
function  on  the  open  interval  (0,a).  Therefore,  0  does  not  converge 
uniformly  to  /kt)  on  the  open  interval  (0,a). 

In  Figure  28  the  phase  of  the  fluctuating  temperature,  O,  is  plotted  for 
0<r<a,  0  =  7r/4,7t/2,3jt/4,7t,  and  Q  =  1,10,100,1000.  The  phase  of  0  is 
not  shown  for  0  =  0  since  0  s  0  for  this  case.  The  phase  of  the  fluctuating 
temperature  field  obtained  using  the  two-term  expansion  for  is  also 
included  for  comparison.  Note  that  in  Eq.  (174)  the  difference  between  the 
solutions  corresponding  to  the  two-term  and  three-term  expansions  for  is 
the  sin(30/2)  term.  Thus,  the  temperature  solution  corresponding  to  the  two- 
term  expansion  for  may  be  separated  in  terms  of  the  0-coordinate  and  the 
radial  coordinate,  i.e.  0  corresponding  to  the  two-term  solution  has  the  form 


Phase  of  the  Fluctuating  Temperature,  o  Phase  of  the  Fluctuating  Terrperature,  o 


100 


Normalized  Radial  Coordinate,  r/a  Normalized  Radial  Coordinate,  xta 


(a)  (b) 


0.00  0.25  0.50  0.75  1.00  0.00  0.25  0.50  0.75  1.00 

Normalized  Radial  Coordinate,  r/a  Normalized  Radial  Coordinate,  r/a 

(c)  (d) 


Figure  28.  Phase  of  the  fluctuating  temperature  as  a  function  of  the 
normalized  radial  coordinate  for  0  =  7r/4,7i/2,37i/4,7i.  (a)  Q  =  1;  (b)  Q=  10; 
(c)  Q  =  100;  (d)  Q  =  1000.  In  the  case  of  the  two-term  expansion,  O  is 
independent  of  0.  Mode  II. 
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0  ~  sin(9/2)x(some  complex-valued  flmction  of  r).  It  follows  that  the  phase 
of  the  fluctuating  temperature  field  corresponding  to  the  two-term  expansion 
for  is  independent  of  0.  This  result  is  evidenced  in  Figure  28. 
Furthermore,  from  Figure  27  it  is  clear  that  the  sin(30/2)  term  in  Eq.  (174) 
makes  a  significant  contribution  to  ©  only  for  r  >  a/2 .  Therefore,  the  phase 
of  the  temperature  solution  corresponding  to  the  three-term  expansion  for 
is  virtually  independent  of  0  for  r<a/2.  For  Q«l,  the  problem  is 
essentially  isothermal,  and  we  observed  that  O  is  nearly  n/l  (plots  not 
included).  Conversely,  for  Q  =  1000,  away  from  the  crack  tip,  the  problem  is 
essentially  adiabatic,  and  O  is  nearly  zero.  At  the  crack  tip,  for  high 
frequencies,  Eq.  (180)  predicts  a  phase  0  =  7t/4,  which  is  evidenced  in 
Figure  28c,  and  28d.  In  summary,  as  Q  increases  from  zero  to  infinity,  the 
phase  decreases  from  7t/2  to  7t/4  at  the  crack  tip,  and  from  7t/2  to  zero  away 
from  the  crack  tip. 

4.3.2  Work  Lost 

The  conversion  of  external  work  into  heat  at  a  point  is  now  examined.  A 
normalized  work  lost  v(/  is  defined  as  vif  =  Aw/(T^K^7^a^/c)  for  Mode  II.  A 
"physical  feel"  for  the  denominator  may  be  developed  as  follows.  A  rod  of 
one  unit  volume  at  an  equilibrium  temperature  7^  is  subjected  to  a  uniaxial 
tension  of  magnitude  x  under  adiabatic  conditions;  then  V  =  —'iTja./C.  While 
holding  T  fixed,  the  rod  is  allowed  to  absorb  heat  from  a  reservoir  at  until  it 
returns  to  its  equilibrium  temperature  2^.  Next,  the  rod  is  unloaded 
adiabatically;  then  V  =  +x7^a/C.  Once  again,  the  rod  is  allowed  to  give  up 
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heat  to  a  reservoir  at  until  the  rod  temperature  becomes  2^ .  The  total  work 
lost  during  this  cycle  is  x^Tja? /C.  By  way  of  example,  for  steel, 

a  =  12xl0”*ir“’  and  C=  3.8x10*  A".  If  we  choose  t  =  100 MPa, 

=  300 J^,  and  k  =  1  (plane  stress),  then  Avi'=  IHJ/tw^  .  Interestingly,  the 
temperature  increase  of  the  rod  due  to  this  conversion  of  work  into  internal 
heat  energy  is  only  AT  =  3x1  O'*  .S'.  For  use  in  the  next  section,  the  total 
entropy  produced  during  this  cycle  is  x^a^/C.  For  the  example  under 
consideration,  As  =  Q.2ZJln^K.  Substituting  Eq.  (162)  into  Eq.  (24)  and 
noting  the  difference  in  notation,  yields 


v|;(r,0)  = -27t 


sin(0/2) 


2  V  r  4Ua  I  16^^“^  ^  A2aJ 


lm(©(r,0)) 


for  0<r  <a,  O<0<7t.  (181) 
where  ©  is  given  by  Eq.  (174).  As  expected,  \|;  =  0  on  0  =  0  because  =  0 


on  0  =  0. 

The  crack-tip  work-lost  field  is  now  examined.  For  a  finite  Q,  expanding 
Eq.  (181)  in  (r/a)  about  r  =  0  and  retaining  upto  the  linear  term,  yields 


vi;  =  jt-\/2(lm  sin^  (0/2) + 


+  71 


-Q  sin^  (0/2)  ''/2  (imfi)  sin(0/2)  sin(30/2) 


r«l  (182) 

where,  as  before,  the  complex  function  ^(q)  is  given  by  Eq.  (178),  and  the 
new  complex  function  B(q)  is  given  by 
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5(a)  =  -^[3-(/Q  +  3(l  +  /Q))e-'®]  (183) 

The  first  term  in  Eq.  (182)  contains  contributions  from  the  first  two  terms  of 
(5^^  in  Eq.  (162),  but  not  from  any  higher-order  terms.  Next,  the  three  terms 
inside  the  square  brackets  in  Eq.  (182)  are  examined:  The  first  term, 
-f2sin^(0/2),  corresponds  only  to  the  first  term  of  in  Eq.  (162);  the 
second  term,  — (3/2-v/2)(lm.(4)sin^(0/2),  contains  contributions  only  from  the 
first  two  terms  of  <3^  in  Eq.  (162);  the  third  term, 
V2(lm5)sin(0/2)sin(30/2),  contains  contributions  only  from  the  first  and 
third  terms  of  in  Eq.  (162).  More  importantly,  if  the  series  expansion 
Eq.  (162)  had  been  carried  out  to  infinity,  none  of  these  three  terms  would 
have  contained  contributions  from  any  of  the  higher-order  terms  in  Eq.  (162). 
Finally,  the  complex  function  B(q)  corresponds  only  to  the  third  term  of 
in  Eq.  (162). 

Next,  the  limit  of  the  work  lost,  \|/,  as  r  0  is  examined;  from  Eq.  (182) 
vj;  =  7l^/2(Im.4)sin^(0/2)  (184) 

Note  that  analogous  to  the  crack-tip  stress  field  in  Eq.  (162)  and  the  crack-tip 
temperature  field  in  Eq.  (177),  the  limit  of  also  depends  upon  the  direction  0 
along  which  the  limit  is  taken.  For  0  =  ti,  ij/  at  the  crack  tip  is  plotted  against 
Q  in  Figure  29.  For  Q  <  10,  v|/  is  plotted  as  a  dashed  line  to  remind  the  reader 
that  these  values  may  be  inaccurate.  The  slope  changes  from  one  as  Q  -»  0  to 
1/2  as  Q  00 . 
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Figure  29.  Normalized  work  lost  at  the  crack  tip,  v);  at  0  =  7t,  as  a  function  of 
the  normalized  frequency  The  dashed  line  indicates  the  range  of  Q  in 
which  the  solution  may  be  inaccurate.  The  slope  decreases  from  one  at  low 
frequencies  to  1/2  at  high  frequencies.  Mode  II. 
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Next,  the  high-frequency  behavior  of  the  work  lost  at  the  crack  tip  is 
examined.  As  Q ^  oo,  Im^(Q) VQ ,  ImB(Q) ->  0,  and  Eq.  (182)  reduces 
to 

M/  =  7i(V2Q-Q(r/fl))  sin^  (e/2) .  r  «l« a  (185) 

or,  in  dimensional  form 


Aw  = 


•yfl  I —  or 

—  —  Vm - 

k 


Kl  sm^{Q/2)  ,  r«l«a  (186) 


Note  that  Aw  ~  ^[o  at  the  crack  tip.  As  in  the  case  of  the  temperature  V  in 
Eq.  (180),  Aw  in  Eq.  (186)  corresponds  only  to  the  square-root  singular  stress 
field  in  Eq.  (162).  Therefore,  for  a  fixed  Kj,  and  o,  Eq.  (186)  describes  the 
universal  high-frequency  work  lost  field  at  the  crack  tip  for  all  Mode  II 
problems. 

Returning  to  Eq.  (181),  in  Figure  30  vj/  is  plotted  for  0<r<a, 
0  =  0,7t/4,7t/2,3Tc/4,7t,  and  0  =  1,10,100,1000.  Once  again,  even  though  it 
may  be  quite  inaccurate,  v|;  is  plotted  for  O  =  1  to  give  a  "feel"  for  it.  Note  that 
the  Figure  for  O  =  1000  has  a  different  scale.  The  work  lost  obtained  using 
the  two-term  expansion  for  is  also  included  for  comparison.  The  results 
obtained  with  the  two-term  and  three-term  expansions  are  virtually 
indistinguishable  for  all  four  values  of  Q.  This  result  provides  some 
justification  for  our  assumption  to  truncate  the  Laurent  series,  Eq.  (162),  at 
three  terms.  The  work  lost  is  essentially  confined  to  the  interior  of  the  disk 
r  =  a.  This  provides  an  a  posteriori  justification  for  replacing  the  actual 


Normalized  Work  Lost,  \|/  Normalized  Work  Lost, 
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(C)  (d) 

Figure  30.  Normalized  work  lost  as  a  flmction  of  the  normalized  radial 
coordinate  for  0  =  0, 7i/4,7i/2,37t/4,7t.  (a)  Q  =  l;  (b)  fi  =  10;  (c)  £2  =  100; 
(d)  £2  =  1000.  The  work  lost  corresponding  to  the  two-term  and  the  three-term 
expansions  is  virtually  indistinguishable.  The  work  lost  is  essentially 
confined  to  the  interior  of  the  disk  r  =  a.  Mode  11. 
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g(r,Q)  -  which  is  very  small  any  way  -  with  zero  for  r  >  a.  Asa  check  on  our 
solution,  we  note  that  in  the  exact  solution  on  the  plane  x  =  0,  V  =  0,  and 
therefore  \\f  =  0.  In  the  approximate  solution  the  only  point  that  coincides 
with  the  plane  x  =  0  is  (r  =  a,0  =  ti)  .  In  Figure  30  at  (r  =  a,0  =  7t) ,  ij/  is  zero 
for  all  Q,  thus  giving  credence  to  our  approximate  solution.  In  physical  terms, 
at  low  frequencies  essentially  isothermal  conditions  are  obtained,  and  the 
work  lost  is  nearly  zero  throughout.  This  case  is  evidenced  by  Q  =  1  in 
Figure  30.  Conversely,  at  high  frequencies,  away  from  the  crack  tip, 
essentially  adiabatic  conditions  are  obtained,  and  again  the  work  lost  is  zero. 
Explicitly,  from  Eq.  (181),  for  any  r  e(0,a),  limv(;(r,0)  =  0.  From 

Q~>oo 

Eq.  (186)  we  see  that  \|/  ~  V©  at  the  crack  tip.  Therefore,  as  in  the  case  of  the 
temperature,  v|;  does  not  converge  uniformly  to  its  adiabatic  limit,  x|/(r,0)  =  0, 
on  the  open  interval  (0,a). 


4.3.3  Entropy  Produced 

The  production  of  entropy  at  a  point  is  now  examined.  A  normalized 
entropy  produced  is  defined  as  x  =  As/(T^K^a^/c).  A  physical  feel  for  the 
denominator  was  discussed  in  the  preceding  section.  Rederiving  Eq.  (27)  in 
polar  coordinates  and  noting  the  difference  in  notation,  gives 


a©  a© 

a(r/a)  d{r/a) 


1  a©  a© 

{r/af  'mle 


(187) 


where  ©  is  given  by  Eq.  (174). 
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First,  the  crack-tip  entropy  field  is  examined.  For  a  finite  Q  (finite  /), 
expanding  Eq.  (187)  about  r  =  0  and  retaining  only  the  first  term  yields 

nAA  a 


%  =  ■ 


r«l 


(188) 


4Q  r’ 

where  .(4(q)  is  given  by  Eq.  (178).  As  in  the  case  of  ©  in  Eq.  (177),  %  in 
Eq.  (188)  contains  contributions  from  the  first  two  terms  of  a**  in  Eq.  (162), 
but  not  from  any  higher-order  terms.  Observe  that  %  exhibits  a  1/r  type 
singularity  at  the  crack  tip,  and  is  independent  of  0.  A  normalized  entropy 
intensity  factor  is  defined  as 


=H --x1= 


AA 

2Q. 


(189) 

r-^°yna’'J  in 

Figure  31  shows  the  dependence  of  Sjj  on  For  Q  <  10,  Sj,  is  plotted  as  a 
dashed  line  to  indicate  that  these  values  may  be  inaccmate. 

Next,  the  high-frequency  behavior  of  the  crack-tip  entropy  field  is 
considered.  For  O  oo ,  A  A  /Q  ->  2  so  that  therefore 


Tz  a 


2  r 


r«l«a 


(190) 


or  in  dimensional  form 


.  tikV  KI  . 

As’  = - —  ,  r«l«a  (191) 

ClTzr 

In  the  spirit  of  LEFM,  a  dimensiorial  enixoTpy  intensity  factor  is  defined  as 

2  2 

M,j  =  lim(27ir  Ay)  =  ^  j  ^2) 


NormalLzied  Entropy  Intensity  Factor,  S 


no 


so  that 

M 

As  =  —^,  r«l«a  (193) 

Inr 

Note  that  in  the  high-frequency  limit,  x  in  Eq.  (190)  is  frequency  independent 
(cf.  the  temperature  0  ~  Vq  ,  Eq.  (179),  and  the  lost  work  i)/  ~ 

Eq.  (185)).  As  in  the  case  of  the  temperature  V  in  Eq.  (180)  and  the  work  lost 
Aw  in  Eq.  (186),  As  in  Eq.  (191)  corresponds  only  to  the  square-root  singular 
stress  field  in  Eq.  (162).  Therefore,  for  a  fixed  Ku,  Eq.  (191)  describes  the 
universal  high-frequency  crack-tip  entropy  field  for  all  Mode  II  problems. 

In  Figure 32,  x  plotted  for  0<r<a,  0  =  O,7i/4,7t/2,37T/4,7r,  and 
Q  =  1,10,100,1000.  The  entropy  produced  using  the  two-term  expansion  for 
is  also  included  for  comparison.  Once  again  the  results  obtained  with  the 
two-term  and  three-term  expansions  are  virtually  indistinguishable  for  all  the 
Q  shown,  lending  further  credence  to  our  decision  to  truncate  the  series  at 
three  terms.  Note  that  the  normalization  of  Aw  and  As  was  done  in  such  a 
manner  that  setting  v|;  =  x  corresponds  to  Aw  =  7^As,  i.e.  one  unit  of 
normalized  entropy  produced  corresponds  to  one  unit  of  normalized  work  lost. 
Accordingly,  \}/  (Figure  30)  and  x  (Figure  32)  are  plotted  using  identical  scales 
(except  Figure  30d).  As  in  the  case  of  xj/,  as  r  ->  a,  x  is  negligibly  small  for 
all  Q.  This  provides  further  justification  for  replacing  g(r,0)  -  which  was 
very  small  any  way  -  with  zero  fox  r>  a.  It  is  instructive  to  compare  the  lost 
work  (Figure  30)  and  the  entropy  produced  (Figure  32).  The  lost  work  is 
identically  zero  on  0  =  0  (in  front  of  the  crack),  increases  monotonically 


Normalized  Entropy  Produced,  Normalized  Entropy  Produced,  ^ 
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Normalized  Radial  Coordinate,  r/a  Normalized  Radial  Coordinate,  xta 


(a)  (b) 


Normalized  Radial  Coordinate,  ria  Normalized  Radial  Coordinate,  r/a 

(c)  (d) 


Figure  32.  Normalized  entropy  produced  as  a  function  of  the  normalized 
radial  coordinate  for  0  =  0, 7t/4,7i/2,37t/4,7t.  (a)  Q  =  l;  (b)  Q  =  10;  (c) 

Q  =  100;  (d)  Q  =  1000.  The  entropy  produced  goes  to  infinity  as  (l/r)  at  the 
crack  tip.  The  entropy  produced  corresponding  to  the  two-term  and  the  three- 
term  expansions  is  virtually  indistinguishable.  Mode  II. 
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with  0,  and  becomes  maximum  on  0  =  7t  (crack  faces).  Close  to  the  crack-tip, 
the  exact  opposite  is  true  for  x-  (This  does  not  come  as  a  surprise.  In 
Section  3.1.2  it  was  shown  that  in  the  case  of  an  Euler-Bemoulli  beam,  x  =  0 
on  the  adiabatic  boundaries  and  maximum  on  the  neutral  axis,  whereas  v|/  s  0 
on  the  neutral  axis  and  maximum  on  the  boundaries).  In  Figure  32  for 
Q  =  100  and  Q  =  1000  there  is  a  fairly  sharp  cross-over  visible  at  rja  =  0.2 
and  r/a  =  0.08,  respectively.  For  rja  greater  than  the  cross-over  point,  %  is 
lowest  on  0  =  0  and  increases  Avith  0.  For  any  r  €  (0,n),  in  the  high  frequency 
limit,  Eq.  (187)  reduces  to  limx(r,0)  =  O.  On  the  other  hand,  from 

Q,-^co 

Eq.  (190),  X  is  unbounded  at  the  crack  tip.  Therefore,  as  in  the  case  of©  and 
X  does  not  converge  uniformly  to  the  adiabatic  limit  on  the  open  interval 
(0,a). 

4.3.4  Total  Work  Lost 

The  total  work  lost  throughout  the  plate  (of  unit  thickness)  per  cycle  is 
given  by 

AW  =  f  f  Awdxcfy  (194) 

J—tO  J— 00 

However,  in  the  previous  section  it  was  argued  that  virtually  all  of  the  external 
work  is  lost  in  the  interior  of  the  disk  r<a.  Accordingly,  Eq.  (194)  is 
approximated  with 

AW=2j  iBj^Awrdr  (195) 

where  the  factor  of  2  accounts  for  the  two  crack  tips.  With  reference  to 
Figure  30,  it  is  clear  that  the  solution  for  Aw  obtained  using  the  two-term 


113 


expansion  for  is  virtually  identical  to  the  solution  obtained  using  the 
three-term  expansion.  Therefore,  for  simplicity  the  two-term  solution  will  be 
used.  (This  solution  is  obtained  by  omitting  the  sin(30/2)  terms  in  Es.  (174) 
and  (181).)  We  now  define  a  normalized  total  work  lost,  T  =  AW/W,  where 
W  =  I c).  (Note  that  the  term  in  parenthesis  is  the  term  that 

arose  naturally  in  the  discussion  of  the  work  lost,  \j;.  Once  again,  the  factor  of 
2  accounts  for  the  two  crack  tips.)  Then 

^  =  2%  Im{^  -  V/O  +  e~'^  -  3)(^  -  ^  cosh( ^  sinh(-v/^)) j} 

(196) 

In  Figure  33,  T'  is  plotted  as  a  function  of  Q.  For  Q  <  1 0  we  have  plotted 
'P  as  a  dashed  line  to  indicate  that  these  values  may  be  inaccurate.  As 
expected,  as  Q^O,  essentially  isothermal  conditions  are  obtained  and 
T  ->  0.  As  Q  ->•  00,  away  from  the  crack  tip,  essentially  adiabatic  conditions 
are  obtained,  and  once  again  ^->-0.  The  maximum  work  lost  occurs  at 
approximately  Q  =  10. 

The  exact  solution  for  the  total  work  lost  in  an  Euler-Bemoulli  beam  of  a 
rectangular  cross-section  may  be  obtained  by  substituting  Eq.  (41)  into 
Eq.  (3),  yielding 


W 


192  A  Q 

odd 


(beam) 


(197) 
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Figure  33.  Normalized  total  work  lost  as  a  fimction  of  the  normalized 
frequency.  The  dashed  line  indicates  the  range  where  the  solution  may  be 
inaccurate.  Mode  II.  The  corresponding  results  for  an  Euler-Bemoulli  beam 
are  also  included  for  comparison  (dotted  line).  The  close  resemblance 
between  the  two  elastothermodynamic  responses  comes  as  a  pleasant  surprise. 
The  peak  occurs  at  virtually  the  same  normalized  frequency,  Q  =  10. 
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where  W  =  {‘KEa^T^/c)w^,  is  the  maximum  stored  elastic  energy  during 

the  cycle,  E  is  the  Young's  modulus,  Q  =  GiCh^/k  (cf.  Q  =  (i>Ca^/k  for  the 
Griffith  crack),  and  h  is  the  total  height  of  the  beam  in  the  plane  of  the 
vibrations.  This  result  was  originally  derived  by  Zener  (1937)  using  a 
different  approach.  Moreover,  Eq.  (197)  has  been  repeatedly  verified  by 
several  experimentalists,  e.g..  Berry  (1955)  or  Nowick  and  Berry  (1972, 
Ch.  17).  The  exact  solution  for  the  beam  (dotted  line)  and  our  approximate 
solution  for  the  crack  are  compared  in  Figure  33.  In  spite  of  the  obvious 
differences  in  the  essential  nature  of  the  two  boimdaiy  value  problems,  there  is 
a  close  resemblance  between  their  elastothermodynamic  responses.  The  peak 
occurs  at  virtually  the  same  normalized  frequency  Q;  the  semi  crack  length  (a) 
and  the  beam  height  (h)  play  the  same  role  in  the  definition  of  Q. 
Furthermore,  there  is  a  similarly  close  resemblance  between  the  approximate 
elastothermodynamic  response  of  the  Griffith  crack  and  the  exact  response  of 
a  periodic  array  of  slabs  (Section  3.4.1),  the  symmetric  three-layer  plate 
(Section  3.4.2),  and  the  composite  sphere  (Section  3.4.3).  This  comparison 
lends  further  credence  to  the  simplifying  assumptions  introduced  in  obtaining 
the  approximate  solution  to  the  Mode  II  crack  problem. 


116 


4.4  MODE  I 

4.4.1  Temperature  Field 

A  Mode  I  Griffith  crack  subjected  to  a  biaxial  tension  is  shown  in 
Figure  34a.  The  hydrostatic  stress  field  is  given  by  Eq.  (156),  and  is  an  even 
function  of  x  and  y.  The  temperature  field  is  governed  by  Eqs.  (159) 
and  (160).  Since  the  operator  in  Eq.  (159)  is  symmetric  in  x  and  y,  F  is  an 
even  function  of  x  and  y.  Moreover,  the  crack  faces  are  assumed  to  be 
adiabatic.  Therefore,  the  entire  planes  x  =  0  and  y  =  0  are  adiabatic.  Thus,  it 
is  only  necessary  to  consider  the  positive  quarter  plane  (x  >  0,y  >  O)  with  the 
adiabatic  boundary  conditions 


dV 

(1)  —  =  0  x  =  0,0<y<oo 
dx 

dV 

(2)  —  =  0  y  =  0,0<x<oo 
dy 


(198) 


This  is  shown  in  Figure  34b.  Attempts  to  solve  this  heat  conduction  problem 
also  led  to  intractable  mathematics  so  an  approximate  solution  is  again  sought. 

First,  an  approximation  to  is  considered.  As  in  the  case  of  Mode  II, 
we  introduce  a  polar  coordinate  system  (r,0)  at  the  right  crack  tip  z  =  a  (see 
Figure  23),  and  expand  Eq.  (156)  in  a  Laurent  series  in  z/2a  about  the  point 
z  =  a.  Retaining  the  first  three  terms  we  obtain 


=  2kcs\ 


cos(0/2) 


-lcos(3e/2)f:f 

16  V  2u 


,3/2 


0<r  <2a 


(199) 
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Figure  34.  (a)  A  Griffith  crack  of  length  la  in  Mode  I  loading  with  a 
rectangular  coordinate  system  (jc,y).  The  crack  faces  are  assumed  to  be 
adiabatic.  Because  is  an  even  function  of  x  and  y,  the  planes  x  =  0  and 
y  =  0  are  adiabatic,  (b)  The  reduced  quarter-plane  heat-conduction  problem. 
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In  Figure  35  the  exact  contours  (using  Eq.  (156))  are  plotted  for  the  region 
0<x<2a,  0<y<a,  and  in  Figure 36  the  exact  and  its  three-term 
approximation,  Eq. (199),  are  compared  for  0<r/a<2  (cf.  0<r/a<l  for 
Mode  II).  The  comparison  is  quite  good.  A  similar  comparison  was  carried 
our  for  a  four-term  expansion;  it  was  found  that  the  fourth  term  does  not  yield 
any  significant  improvement.  Therefore,  the  series  was  truncated  at  three 
terms.  Note  that  in  Mode  I,  the  exact  =  0  on  the  crack  faces  and,  for  any  r 
increases  to  its  maximum  value  on  0  =  0.  The  exact  opposite  is  true  for 
Mode  II:  the  exact  is  maximum  on  the  crack  faces  and,  for  any  r, 
decreases  to  zero  on  0  =  0  (cf  Figure  25).  As  a  visual  aid,  one  might  say  that 
a  Griffith  crack  tends  to  localize  the  hydrostatic  stress  field  outside  the  crack 
tips  in  Mode  I,  and  inside  the  crack  tips  in  Mode  II.  Returning  to  Figure  36  in 
spite  of  the  singularity  at  the  other  crack  tip  (r  =  2a,  0  =  7t),  the  comparison  is 
considered  to  be  surprisingly  good  out  to  r  <2a.  This  may  be  attributed  to 
the  fact  that  =  0  on  the  crack  faces,  0  <  r  <  2a,  0  =  ±7t . 

The  approximations  introduced  in  the  paragraph  containing  Eq.  (162)  are 
also  used  for  Mode  I.  The  resulting  "associated  heat-conduction  problem"  is 
shown  in  Figure  37.  Next,  we  consider  the  heat  generation  term  ^(/■,0). 
Inside  the  region  0^0 <7t,  0<r<2a,  the  exact  g(r,Q)  (Eqs. (156) 
and  (160))  is  approximated  by  the  three-term  expansion  (Eqs.  (199)and  (160)). 
Now,  as  r  ->  00  (for  any  0),  the  exact  /kct  2 .  As  r-^  2a  (for  any  0), 
the  approximate  is  very  nearly  two  (see  Figure  36). 
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Figure  35.  Mode  I  normalized  hydrostatic  stress  contours,  g^^/kc.  (x/a  =  1 
corresponds  to  the  right  crack  tip.  The  full  field  may  be  constructed  by  using 
the  even  symmetry  of  the  problem.  The  range  of  the  contours  is  0,1  (lower 
left)  to  2.0  (right  of  crack  tip)  with  increments  of  0.1.) 


Normalized  Hydrostatic  Stress,  oj^  /  kct 
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Figure  36.  Comparison  between  the  exact  and  its  three-term  expansion 
for  0  =  0,  7i/4,7:/2,37i/4,7i.  Model. 
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i  gir,<S)=g(2a,Q) 


r=2a 


Figure  37.  The  "associated  heat  conduction  problem"  with  a  nonzero  heat 
generation,  ^(r,0)  =  g(2a,0),  outside  the  disk  of  radius  2a.  Mode  I. 
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A  jump  discontinuity  in  stress  would  lead  to  steep  gradients  in  temperature 
which,  in  turn,  would  lead  to  spurious  entropy  production  Eq.  (11).  In  order 
to  avoid  this  situation,  in  the  associated  heat-conduction  problem,  g(r,Q)  is 
set  equal  to  g(2a,Q)  in  the  exterior  of  the  half-disk,  i.e.  in  O<0<7i:, 
2a<r  <<x>  (cf.  Mode  II  where  we  set  g(r,Q)  sO,  inO<0<7r,a<r<  oo).  As 
in  the  case  of  Mode  II,  once  the  solution  to  the  associated  heat-conduction 
problem  in  y  >  0,  the  solution  in  the  left  quarter-plane,  x  <  0,  y  >  0,  will  be 
discarded;  the  solution  in  the  right  quarter  plane,  x  >  0,  y  >  0,  will  be  taken  as 
an  approximation  to  the  exact  solution  of  the  problem  shown  in  Figure  34b. 
The  solution  in  the  full  plane.  Figure  34a,  can  then  be  constructed  using  the 
even  symmetries  of  F. 

The  associated  heat-conduction  problem  for  Mode  I  is  given  by 

|^+i^+-L^-£toF+ig(r.e)  =  0,  osesn,  0<r<«.  (200) 
dr  r  dr  r  50  k  k 

dV 

—  =  0  at  0  =  0, 7C,  0<r<oo  (201) 

50 

where 


J  -{i(Q)(xT^a^{r,Q)  O<0<7t,  0<r<2a 

|-(/co) <5^^ (2a, 0)  O<0^7r,  2a<r<oo 


(202) 


An  exact  solution  to  this  problem  was  obtained  by  using  an  integral  transform 
method  analogous  to  that  used  for  Mode  II;  the  solution  is 
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0(r.8)  =  -V2in2;^25W. 

n=0 

for  0  <  r  <  2a  (203) 

where  N„  =  n  for  n  =  0  and  7i/2  for  «  =  1,2,3,...,  and  a  normalized 
fluctuating  temperature  ©  for  Mode  I  has  been  defined  as  ©  =  F/(aK2^a/G) 
(cf.  Mode  II,  where  ©  =  fy(TK7^a/C)).  Unfortunately,  the  integrals  could  not 
be  evaluated  in  a  closed  form.  Therefore,  numerical  integration  was  used. 
The  Bessel  functions  and  the  associated  integrals  were  evaluated  using  version 
2.0  of  the  IMSL  subroutines  CBINS,  CBKS,  and  QDAG  with  the  parameters 
ERRREL  set  equal  to  0.005  and  ERRABS  set  equal  to  0.0.  The  improper 
integrals  were  evaluated  by  taking  the  upper  limit  to  be  a  large  constant  in 
place  of  infinity.  It  was  found  that  an  upper  limit  of  10  was  sufficiently 
accurate  for  graphing  purposes  for  the  range  of  frequencies  considered. 
Furthermore,  it  was  found  that  the  first  eleven  terms  in  the  summation  in 
Eq.  (203)  were  sufficient. 

Unlike  Eq.  (174),  Eq.  (203)  could  not  be  reduced  to  a  closed  form. 
Therefore,  an  expression  for  the  crack-tip  temperature  field  could  not  be 
derived  (cf.  Eq.  (177)  for  Mode  II).  However,  it  was  proven  from  Eq.  (203) 
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that  for  any  finite  fi,  the  temperature  is  bounded  at  the  crack  tip.  So  far,  for 
maximum  accuracy,  the  solution  has  been  constructed  in  the  half-disk,  r  <  2a ; 
the  approximation  to  g(r,B)  starts  at  r  =  2a.  However,  from  the  symmetries 
of  ^(x,^)  it  is  obvious  that  only  the  solution  in  the  half-disk  r  <a  is  needed 
(see  Figure  34b).  Accordingly,  in  Figure  38  the  magnitude  of  0  is  plotted  for 
0^r/a<l,  0  =  0,7i/4,7i/2,37t/4,7i,  and  Q  =  1,10,100,1000.  For  the  reasons 
discussed  in  the  paragraph  following  Eq.  (163),  it  is  recognized  that  the 
approximate  solution  may  not  be  accurate  for  Q<10.  Nevertheless,  0  is 
plotted  for  Q  =  1  to  give  a  "feel"  for  it.  The  results  using  a  two-term 
expansion  for  are  also  shown  for  comparison.  For  a  given  r  the 
temperature  magnitude  decreases  monotonically  from  a  maximum  at  0  =  0  to 
a  minimum  at  0  =  tt  (cf.  Figure  27  where  it  is  the  exact  opposite).  It  is  clear 
that  for  Q  =  1,  heat  has  plenty  of  time  to  diffuse,  and  essentially  isothermal 
conditions  are  obtained.  Now  examine  the  other  extreme,  Q  =  1000.  If  truly 
adiabatic  conditions  are  obtained  then  0  =  — Cit  /kct.  By  overlaying  Figure  36 
and  Figure  3  8d,  drawn  to  the  same  scale,  we  observed  that  beyond  some 
distance  away  from  the  crack  tip,  essentially  adiabatic  conditions  are  in  fact 
achieved  at  =  1000.  Due  to  the  limitations  of  evaluating  0  numerically,  we 
were  only  able  to  calculate  0  for  r/a  >0.01. 

In  Figure  39  the  phase  of  the  fluctuating  temperature,  O,  is  plotted  for 
0<r<a,  0  =  0,7i/4,7t/2,37i:/4,7i,  and  0  =  1,10,100,1000.  The  phase  of  the 
fluctuating  temperature  field  obtained  using  the  two-term  expansion  for  is 
also  included  for  comparison.  Note  that  in  contrast  to  Eq.  (174),  in  Eq.  (203) 
the  temperature  solution  corresponding  to  the  two-term  expansion  for 


N^bgnitude  of  the  Normalized  Magnitude  of  the  Normalized 

Fluctuating  Tenqjerature,  ||©i|  Fluctuating  Temperature,  ||G|| 


125 


Normalized  Radial  Coordinate,  r/a  Normalized  Radial  Coordinate,  xia 
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Figure  38.  Magnitude  of  the  normalized  fluctuating  temperature  as  a  function 
of  the  normalized  radial  coordinate  for  0  =  0, 7i:/4,7i/2,37i/4 ,7i .  (a)  Q  =  1;  (b) 
Q  =  10;  (c)  Q  =  100;  (d)  Q  =  1000.  The  fluctuating  temperature  is  boimded  at 
the  crack  tip.  Model. 
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Figure  39.  Phase  of  the  fluctuating  temperature  as  a  function  of  the 
normalized  radial  coordinate  for  0  =  0, 7t/4,7i/2,37i/4,7t.  (a)  Q  =  l;  (b) 

Q  =  10;  (c)  Q  =  100;  (d)  Q  =  1000.  Mode  I. 


127 


can  not  be  separated  in  terms  of  the  9-coordinate  and  the  radial  coordinate. 
There  is  a  much  stronger  9-dependence  in  Mode  I  (Figure  39)  as  compared  to 
Mode  II  (Figure  28).  Observe  that  for  large  frequencies,  Q»1999,  the 
problem  is  essentially  adiabatic  away  from  the  crack  tip,  and  the  fluctuating 
temperature  is  out  phase  with  the  applied  stress  by  approximately  n. 
Conversely,  for  small  frequencies,  O  » 1,  the  problem  is  essentially  isothermal 
and  once  again,  the  fluctuating  temperature  is  out  phase  with  the  applied  stress 
by  approximately  iz.  As  the  frequency  increases,  spatial  variations  in  O 
become  increasingly  localized  near  the  crack  tip. 


4.4.2  Work  Lost . 

For  Model,  a  normalized  work  lost  \\>  is  defined  as 
i|/  =  fc)  (cf.  \\i  =  Aw/(T^K^7^aVc)  for  Mode  II).  Substituting 

Eq.  (199)  into  Eq.  (24)  and  noting  the  difference  in  notation,  yields 


\|/(r,9)  =  -271 


cos(9/2) 


2  \  r  4  V2a 


■icos(39/2)f 
10 


3/2 


Im©(r,9) 


for9<r<2cr.  (294) 

where  0  is  given  by  Eq.  (293).  As  expected,  in  view  of  Eq.  (24),  on  9  =  tc, 
v[f  s  0  because  =  0  (cf.  Mode  II  where  \|;  =  9  on  9  =  0).  Since  the 
temperature  is  boimded  at  the  crack  tip,  it  follows  from  Eq.  (204)  that  the 
work  lost  exhibits  a  l/Vr  singularity  at  the  crack  tip.  A  dimensionless  work- 
lost  intensity  factor  v|/^  is  defined  for  Mode  I  as 


r-»0  '0=0 


(205) 
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To  obtain  this  limit,  we  used  a  rational  function  extrapolation  (Press,  1989,  pp. 
83-86).  Figure  40  shows  the  frequency  dependence  of  \\fj  in  the  range 
10  <  <  1000  (The  frequency  range  was  limited  due  to  numerical  difficulties 

in  evaluating  Eq.  (203).). 

In  Figure 41  v|/  is  plotted  for  0<r<a,  0  =  O,7i/4,7r/2,37i;/4,7t,  and 
Q  =  1,10,100,1000.  The  results  obtained  using  a  two-term  expansion  for 
are  also  shown  for  comparison.  For  a  given  r,  the  maximum  \\i  occurs  for 
0  =  0.  There  are  regions  where  \|/  becomes  negative.  A  similar  result  was 
found  in  Section  3.4.1  where  it  was  noted  that  the  only  constraint  imposed  by 
the  second  law  of  thermodynamics  is  that  the  total  entropy  produced 
throughout  the  cracked  plate  (a  closed  thermodynamic  system  which  does  not 
exchange  heat  with  its  surroundings)  during  one  cycle  of  loading  is  either  zero 
or  positive,  i.e.  AS  >  0.  It  follows  from  Eq.  (35)  that  2^ AS  =  AW >  0,  i.e.,  the 
total  work  lost  throughout  the  cracked  plate  during  one  cycle  should  be  zero 
or  positive.  It  then  follows  that  >  0.  From  a  visual  examination  of 

Figure  41,  it  is  apparant  that  this  second-law  constraint  is  clearly  satisfied. 

4.4.3  Entropy  Produced 

A  normalized  entropy  produced  per  unit  volume,  %,  for  Mode  I  is  defined 
as  x  =  Asj{cs\^a^ Ic)  (cf.  %  =  Asl{x\^a^  jc)  in  Mode  II).  Then  %  is  still 
given  by  Eq.  (187),  and  for  brevity  the  expression  resulting  from  the 
substitution  of  Eq.  (203)  into  Eq.  (187)  is  not  given  herein.  It  was  proven  that 
the  entropy  produced  is  bounded  at  the  crack  tip.  In  Figure  42,  %  is  plotted  for 


129 


logjo^ 

Figure  40.  Normalized  work  lost  intensity  factor,  as  a  fimction  of  the 
normalized  frequency  Mode  1. 
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Figure  41.  Normalized  work  lost  \j/  as  a  function  of  the  normalized  radial 
coordinate  for  0  =  0, 7t/4,jt/2,37i/4,7t.  (a)  Q  =  l;  (b)  Q  =  10;  (c)  0  =  100; 
(d)  0  =  1000.  The  work  lost  goes  to  infinity  as  (l/Vr)  at  the  crack  tip. 
Mode  1. 


Normalized  Entropy  Produced,  5^  Normalized  Bitropy  Produced, 


131 


Normalized  Radial  Coordinate,  r/a 

(C) 


Normalized  Radial  Coordinate,  xia 

(d) 


Figure  42.  Normalized  entropy  produced  as  a  function  of  the  normalized 
radial  coordinate  for  0  =  0, 7t/4,7r/2,37i/4,7i.  (a)  Q  =  l;  (b)  Q  =  10;  (c) 

Q  =  100;  (d)  Q  =  1000.  The  entropy  produced  is  bounded  at  the  crack  tip. 
Mode  1. 
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0<r<a,  0  =  O,7i/4,7i/2,37i/4,7i,  and  Q  =  1,10,100,1000.  Note  that  the 
figure  for  Q  =  1000  is  plotted  on  a  different  scale.  The  entropy  obtained  using 
the  two-term  expansion  for  is  also  included  for  comparison.  As  was  done 
in  Mode  II,  vj/  (Figure  41)  and  x  (Figure  42)  are  plotted  using  identical  scales 
(except  Figure  42d).  Sufficiently  close  to  the  crack  tip,  for  a  given  radial 
location,  the  entropy  produced  is  minimum  on  0  =  0  and  increases  to  a 
maximum  on  0  =  tt.  As  a  check  on  the  approximate  solution,  note  that  in  the 
exact  solution  for  Model,  at  the  point  (x  =  0,y  =  0),  dVjdx^Q  and 
dVldy  =  0,  and  therefore  x-0.  In  the  approximate  solution  the  point 
(x  =  0,y  =  0)  corresponds  to  the  point  (r=a,0  =  7c).  In  Figure  42  at 
(/*  =  a,0  =  7i),  X  is  clearly  not  zero  for  Q  =  l,  nearly  zero  for  Q  =  10,  and 
virtually  zero  for  Q  =  100  and  Q  =  1000. 


5.  SUMMARY 


Using  the  linear  one-way  coupled  theory  of  elastothermodynamic 
relaxation,  the  elastothermodynamic  damping  of  composite  and  cracked  media 
was  examined.  Two  descriptions  of  elastothermodynamic  damping  were 
established:  (1)  the  mechanical  description  and  (2)  the  entropic  description. 
An  integral-transform  technique  was  developed  to  solve  for  the  temperatoe 
field  of  a  general  composite  medium  with  arbitrary  heat  generation.  With  this 
solution,  a  general  expression  for  the  damping  of  the  composite  material  was 
derived.  This  general  solution  for  the  elastothermodynamic  damping  was  then 
specialized  to  a  solution  for  a  composite  material  consisting  of  A^-isotropic 
layers  in  a  rectangular,  cylindrical,  and  spherical  coordinate  system  (i.e.  an  N- 
layer  slab,  cylinder,  and  sphere,  respectively).  With  this  specialized  result  the 
elastothermodynamic  damping  of  laminated,  fiber-reinforced,  and  particulate 
composites  was  examined.  Finally,  an  approximate  analysis  was  given  for  the 
temperature  field,  the  work  lost,  and  the  entropy  produced  in  the  vicinity  of  a 
Griffith  crack  for  Modes  I,  II,  and  III. 
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APPENDIX 


To  illustrate  the  physics  of  elastothermodynamic  damping  in  a  simple 
setting,  consider  the  following  quasistatic  problem.  In  analogy  to  a  beam  in 
bending,  two  identical  isotropic  homogeneous  thermoelastic  material  elements 
of  unit  volume,  labeled  A  and  B,  are  taken  through  a  series  of  five  processes  as 
shown  in  the  schematics  of  Figure  A1  and  listed  below.  A  diathermal 
interface  exists  between  the  two  materials  which  may  be  turned  on  and  off  as 
desired. 

process  1-2:  Material  A  initially  at  a  temperature  is  stressed 
adiabatically  from  to  0^=0.  Material  B,  also  initially  at  a 

temperature  T^,  is  stressed  adiabatically  from 

process  2-3:  Materials  A  and  B  transfer  heat  across  the  diathermal 
interface  at  constant  stress  imtil  an  equilibrium  temperature  is  obtained. 

process  3-4:  Material  A  is  stressed  adiabatically  from  =  ct  to  0^;^  =  0. 
Material  B  is  stressed  adiabatically  from  =  -o  to  =  0. 

process  4-5:  Material  A  and  B  transfer  heat  across  the  diathermal 
interface  at  zero  stress  until  an  equilibrium  temperature  is  obtained. 

Note  that  the  system  A-B  is  adiabatic. 
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'j^A 


;ic  of  the  system  A-B  and  the  five  processes  for 
naterial  B  are  similar.  The  curvatures  in  the  processes 
exaggerated.) 
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SOLUTION  USING  THE  LINEAR  ONE-WAY  COUPLED  THEORY 

The  elastothermodynamic  work  lost  for  this  problem  is  now  calculated 
using  linear  one-way  coupled  theory  of  elastothermodynamic  damping,  which 
for  this  simple  example  is  equivalent  to  the  linearization  of  Eq.  (1), 


T,„=T,-T,^ia,„-c,)  (A3) 

where  the  subscripts  denote  the  process  number.  Using  Eq.  (A3)  and  noting 
that  =Tf  ={t,^  +  T2^)l2,  and  =r/ =(r/  +  r/)/2  (the  superscripts 
denote  the  material),  the  temperatures  of  both  materials  A  and  B  at  the 
completion  of  each  process  were  calculated  and  are  given  in  Table  Al.  Note 
that  the  temperature  of  both  materials  at  the  completion  of  the  stress  cycle 
(12345)  is  the  same  as  the  temperature  at  the  beginning  of  the  cycle,  .  Thus, 
within  the  linear  theory,  the  system  A-B  has  imdergone  a  thermodynamic 
cycle. 
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Table  A1  Process  temperatures  (using  Eq.  (A205)) 


MateriaM 

Material  B 

II 

II 

r,  =  r,-r,acT/c 

II 

II 

T4  =  T„+T^aGlC 

r,  =  r,-7;a<T/c 

II 

E-T 

II 

^4“ 

The  external  work  performed  on  a  material,  JV,  during  the  process 
i-^i  +  1  is  given  by 

=  fads  =  ra(c/£”“''  +ds^^)=  ra—-GadT  (A4) 

where  the  fact  that  the  mechanical  strain  8"'“'’  =  ojE ,  where  E  is  the  Yoimg's 
Modulus,  and  the  thermal  strain  =a(2’^., -Tj),  has  been  used.  Using 
Eq.  (A4)  and  the  temperatures  given  in  Table  Al,  the  external  work  performed 
on  materials  A  and  B  during  each  process  was  calculated  and  is  given  in 
Table  A2.  The  net  work  performed  on  the  system  during  the  stress  cycle 
(12345),  AW,  is  the  sum  of  the  work  done  during  each  process  for  both 
materials.  Using  Table  A2,  this  sum  is  equal  to  APP  =  27^  a^a^/C.  (For  a 
numerical  example,  take  steel  where  a  =  12x1  O'*  .S’”'  and  C  =  3.8x10*  j/m^K. 
If  c  =  100 MPa  and  =  300.^,  then  AW  =  227 J  for  a  unit  volume.)  It  is 
instructive  to  calculate  AW  using  the  first  law  of  thermodynamics,  instead  of 
using  Eq.  (A4)  explicitly.  Taking  both  materials  as  our  thermodynamic 
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system,  the  first  law  may  be  written  AC/  =  AW+AQ.  Here,  AQ  =  0;  since 
=  =  T^  =T^,  and  at  the  completion  of  the  stress  cycle  (12345) 

both  materials  are  imstressed,  it  follows  that  AU  =  0  as  well.  Thus,  AfF  =  0, 
contradicting  the  previous  result.  Therefore,  the  linear  one-way  coupled 
theory  violates  the  first  law  of  thermodynamics  (see  the  note  after  Eq.  (Al)). 
This  violation  is  tolerated  since  the  discrepancy  between  the  actual  and  the 
approximate  is  negligible.  (This  fact  will  be  exhibited  in  the  next 
section.) 

Table  A2  Work  performed  during  each  process  (linear  theory) 


Material 


Material  B 


2£  “  2C 


2E  °  2C 


^23  =  T, 


aW 


^23  =  I 


a  a 


^2  ^2^2 
W  --—  +  T  -  - 

"34 


2E  "  2C 


_2 

2E  2C 


^45=0 


«^45=0 


Within  the  linear  one-way  coupled  theory,  the  work  lost  may  also  be 
calculated  using  the  entropy  produced.  The  entropy  change  of  each  material, 
AS,  during  the  process  z  ^  /  + 1  is  given  by 
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f'+i  dT  T  , 

AS„„=C|,  Y  =  Cln^  (A5) 

Using  the  temperatures  given  in  Table  Al,  these  entropies  were  calculated  and 
are  given  in  Table  A3.  To  first  order  in  the  parameter  aa/C,  the  total  entropy 
production  during  the  stress  cycle  (12345)  is  AiS'  =  2a^a^/C.  Thus, 
AS  =  AWfT^  or  AW=TgAS.  In  Section  3.1  this  result  is  proven  using  the 
linear  one-way  coupled  theory,  for  the  general  problem  of  a  composite 
material  made  up  of  several  subregions,  each  of  which  may  be  anisotropic, 
and  whose  interfaces  may  be  thermally  imperfect,  stressed  time-harmonically. 
Even  though  the  system  is  adiabatic,  the  quantity  T^AiS  is  viewed  as  the  "heat 
increment"  per  stress  cycle  since,  as  the  exact  solution  will  indicate,  in  reality 
the  temperature  of  the  system  at  the  completion  of  the  stress  cycle  increases  as 
though  a  quantity  of  heat  had  been  transferred  to  the  system. 


Table  A3  Entropy  change  during  each  process  (linear  theory) 


Material 

Material  B 

o 

11 

K) 

II 

0 

A523=-Cln(l-aCT/C) 

AS23  =-Cln(l  +  aa/C) 

O 

II 

0 

II 

AS'45  =  -Cln(l+acT/C) 

AS45=-Cln(l-ac/C) 
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EXACT  SOLUTION 


For  this  simple  example,  the  elastotheraiodynamic  work  lost  may  be 
calculated  exactly.  Since  for  this  problem,  the  stress  is  independent  of  the 
temperature,  these  two  variables  may  be  separated  in  Eq.  (1)  yielding 

dTlT  =  -{alC)do  (A6) 


Integrating  gives 


(A7) 

Using  Eq.  (A7)  and  noting  that  T-^  =T^  and 

-  {rf  +  T^)l2,  the  temperatures  of  both  materials  A  and  B  at  the 
completion  of  each  process  were  calculated  and  are  given  in  Table  A4.  Note 
that  the  temperature  of  both  materials  at  the  completion  of  the  stress  cycle 
(12345)  is  not  the  same  as  the  temperature  at  the  beginning,  T^,  as  was  the 
case  using  the  linear  theory.  Thus  the  system  A-B  does  not  undergo  a 
complete  thermodynamic  cycle.  (For  the  numerical  example  using  steel,  the 
temperature  at  the  completion  of  the  stress  cycle  is  =T^  =  (1.00003)2^..) 

Using  Eq.  (A4)  and  the  temperatures  given  in  Table  A4,  the  external  work 
performed  on  materials  A  and  B  during  each  process  was  calculated  and  is 
given  in  Table  A5.  The  net  work  performed  on  the  system  during  the  stress 
cycle  (12345)  is  equal  to  ISW  =  2CT^  sinh^(aCT/C).  Recall  the  work  lost  using 
the  linear  theory  was  AW  =  2T„  a^a^/C  which  is  the  first  order  approximation 
to  the  exact  work  lost.  (For  the  numerical  example  using  steel,  the  percent 
difference  between  these  two  results  is  0.0000034%,  thus  justifying  the  use  of 
the  linear  theory.)  As  in  the  linear  approximation,  since  AQ  =  0,  the  first  law 
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reduces  to  AU  =  AW.  Since  A17  =  (r/  -  rj + -T,)  =  2(t,^ -T„),  it  is 
easily  shown  that  the  first  law  is  satisfied  as  expected. 


Table  A4  Process  temperatures  (exact,  using  Eq.  (A7)) 


Material/! 

Material  B 

II 

II 

T2  =  T„e-^'^ 

Ty  =  T^  cosh(aa/ C) 

T^  =  T^  cosh(acy/C) 

T,  = 

=  7^  cosh(aa/C) 

T^  =  T^  cosh(aa/C) 

Table  AS  Work  performed  during  each  process  (exact  solution) 


Material  A  Material  .5 


=^+j;[(acT  +  C>-’'>^-C]  =^+r,[(-ac7  +  C)«“'<^-C] 


ITu  =  -|^-r,[o«j+c(l-«'”'")]  W,.  =-^-r4-a<T+C(l-e-”'^)] 

ir„=o 


Using  (A5)  and  the  temperatures  given  in  Table  A4,  the  entropy  change 
for  each  material  was  calculated  for  each  process,  and  is  given  in  Table  A6. 
Unlike  in  the  linear  approximation,  the  elastothermod5mamic  work  lost  can 
not  be  related  to  the  entropy  produced  with  the  formula  AW  =Tg AS,  since  the 
total  entropy  change  is  AS  =  4Clncosh(aa/ C)  ^  AWfT^ . 


Table  A6  Entropy  change  during  each  process  (exact  solution) 


Material  A 


AS,2=0 


AS23  =  Clncosh(aa/C)+acr 
A534=0 

AS45  =  Clncosh(acr/C)-aa 


Material  B 


AS12  =0 


AS23  =  Clncosh(acy/C)-aa 
AS34=0 

AS45  =  Clncosh(a(T/C)  +  aG 
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